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Information  Models  for  Design  Tolerancing:  From  Conceptual  to 

the  Detail  Design 


R.  Sudarsan,  U.  Roy,  Y.  Narahari,  Ram  D.  Sriram,  K.W. Lyons,  N.  Pramanik 
Abstract: 

Tolerance  design  is  the  process  of  deriving  a description  of  geometric  tolerance  specifications 
for  a product  from  a given  set  of  desired  properties  of  the  product.  Existing  approaches  to  tolerance 
analysis  and  synthesis  entail  detailed  knowledge  of  geometry  of  assemblies  and  are  mostly  applicable 
during  advanced  stages  of  design,  leading  to  a less  than  optimal  design.  During  the  design  process  of 
assemblies,  both  assembly  structure  and  associated  tolerance  information  evolve  continuously  and 
significant  gains  can  be  achieved  by  effectively  using  this  information  to  influence  the  design  of 
an  assembly.  Any  pro-active  approach  to  the  assembly  or  tolerance  analysis  in  the  early  design 
stages  will  involve  making  decisions  with  incomplete  information  models.  In  order  to  carry  out 
early  tolerance  synthesis  and  analysis  in  the  conceptual  product  design  stage,  we  need  to  devise 
techniques  for  representing  function-behavior-assembly  models  that  will  allow  analysis  and  synthesis 
of  tolerances,  even  with  the  incomplete  data  set. 

A ‘ function ’ (what  the  system  is  for)  is  associated  with  the  transformation  of  an  input  physical 
entity  into  an  output  physical  entity  by  the  system.  The  problem  or  customer’s  need,  initially 
described  by  functional  requirements  on  an  assembly  and  associated  constraints  on  the  functional 
requirements  defines  the  concept  of  an  assembly.  This  specification  of  functional  requirements  and 
constraints  define  a functional  model  for  the  assembly.  Many  researchers  have  studied  functional 
representation  ( function  based  taxonomy  and  ontology),  function  to  form  mapping,  and  behavior 
representation  ( behavior  means  how  the  system/product  works).  In  a recent  paper,  [68],  we  presented 
a strong  need  for  comprehensive  function-assembly-behavior  (FAB)  integrated  model. 

In  this  report,  we  discuss  extension  of  the  ideas  presented  in  [68]  and  explain  the  integration 
of  function,  assembly,  and  behavior  representation  into  a comprehensive  information  model  (FAB 
models).  To  do  this,  we  need  to  develop  appropriate  assembly  models  and  tolerance  models  that 
would  enable  the  designer  to  incrementally  understand  the  build-up  or  propagation  of  tolerances 
(i.e.,  constraints)  and  optimize  the  layout,  features,  or  assembly  realizations.  This  will  ensure  ease 
of  tolerance  delivery.  In  an  earlier  paper,  [53],  a multi-level  approach  called  Design  for  Tolerance 
[DFT]  process  was  proposed  which  enables  tolerancing  to  be  addressed  at  successive  stages  of  design 
in  an  incremental  fashion  We  also  address  the  effective  use  of  the  FAB  and  DFT  model  for  design 
tolerancing,  starting  from  conceptual  stage  of  the  design  and  continuously  evolving  throughout  the 
entire  design  process  to  the  final  detailed  design.  These  models  can  eventually  lead  to  tolerance  and 
assembly  standards. 
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1 


Introduction 


Tolerancing  is  a critical  issue  in  the  design  of  electro-mechanical  assemblies.  Tolerancing 
includes  both  tolerance  synthesis  and  tolerance  analysis.  There  are  several  research  thrusts 
that  have  been  reported  in  individual  areas  of  assembly  modeling,  tolerance  analysis,  and 
tolerance  synthesis.  But  an  integrated  and  complete  information  model  that  supports  both 
assembly  activities  (i.e.,  assembly  plan  generation)  and  tolerancing  synthesis  and  analysis 
process  uniformly  throughout  the  entire  product  design  process  does  not  yet  exist.  As 
outlined  in  [53],  an  ideal  assembly  model  for  tolerancing  should:  1)  be  closely  coupled  with 
the  design  process;  2)  be  mutable  through  successive  design  stages;  and  3)  be  capable  of 
representing  assembly  and  tolerance  information  at  any  level  of  abstraction. 

The  other  important  attributes  that  have  been  identified  for  an  appropriate  assembly 
model  are  as  follows:  1)  capturing  design  intent;  2)  embedding  different  views  (relational 
view,  location  logic,  etc.)  in  the  assembly  information  model;  3)  enabling  all  assembly 
information  including  tolerancing  to  be  captured  in  the  model;  and  4)  representing  effective, 
integrated  assembly  information  models  throughout  the  design  process. 

Tolerancing  decisions  can  profoundly  impact  the  quality  and  cost  of  electro-mechanical 
assemblies.  Existing  approaches  to  tolerance  analysis  and  synthesis  in  design  entail  de- 
tailed knowledge  of  geometry  of  assemblies  and  are  mostly  applicable  during  advanced  de- 
sign stages,  leading  to  a less  than  optimal  design  process.  During  the  design  process  of 
assemblies,  both  assembly  structure  and  associated  tolerance  information  evolve  continu- 
ously. Therefore,  significant  gains  can  be  achieved  by  effectively  using  this  information  to 
influence  the  design  of  an  assembly.  Motivated  by  this,  we  identified  and  explored  two  goals 
for  future  research  that  we  believe  can  enhance  the  scope  of  tolerancing  for  the  entire  design 
process  [53,  54].  The  first  goal  addresses  advancing  tolerance  related  decisions  to  the  earliest 
possible  design  stages.  This  issue  raises  the  need  for  effective  representation  of  tolerancing 
information  during  different  design  stages  and  for  effective  assembly  modeling.  The  second 
goal  addresses  appropriate,  synergistic  use  of  available  methods  and  best  practices  for  toler- 
ance analysis  and  synthesis,  at  successive  design  stages.  Pursuit  of  these  goals  lead  us  to  the 
definition  of  a multi-level  approach  that  enables  tolerancing  to  be  addressed  at  successive 
design  stages  in  an  incremental  fashion.  The  resulting  design  process  called  the  design  for 
tolerance  [53]  process  integrates  three  important  domains:  1)  design  activities  at  successive 
design  stages;  2)  assembly  models  that  evolve  continuously  through  the  design  process;  and 
3)  methods  and  best  practices  for  tolerance  analysis  and  synthesis. 

Product  functionality  is  another  important  factor  to  be  considered.  It  is  clear  that 
any  intelligent  decision  during  the  design  of  product  modeling  cannot  be  made  without  the 
knowledge  of  product  functions.  The  information  model  for  assembly  must  include  functional 
and  behavioral  characteristics  of  its  component  parts. 

Setting  tolerances  for  parts  of  an  assembly  has  always  been  a critical  design  decision.  The 
design  of  tolerances  for  component  parts  and  assembly  should  be  cost  effective  and  adequate 
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to  ensure  the  required  performance  specifications  for  an  assembly.  An  inappropriate  choice 
of  tolerances  can  result  in  low  quality  products,  extensive  or  difficult  manufacturing  steps, 
or  both. 

For  complex  designs,  tradition,  trial  and  error,  or  intuition  frequently  determine  toler- 
ances. A common  method  employed  by  designers  is  to  select  the  dimensions  of  a part  that 
are  considered  important,  and  then  specify  tightest  tolerances  that  manufacturing  processes 
can  achieve.  This  unnecessarily  overburdens  the  manufacturing  facilities  without  ensuring 
optimality  of  the  design.  We  need  to  design  tolerances  that  help  in  the  rational  choice  based 
on  considerations  of  cost,  sensitivity,  and  performance  specifications.  However,  we  should 
remember  that  current  assembly  analysis  tools  or  tolerance  analysis/synthesis  tools  have 
always  been  utilized  in  the  detailed  design  phase  of  the  product  development  process. 

In  this  report  our  main  aim  is  to  develop  an  integrated  comprehensive  function,  assembly 
(artifact)  and  behavior  (FAB)1  model.  To  do  this,  we  need  to  develop  appropriate  assembly 
models  and  tolerance  models  that  would  enable  the  designer  to  incrementally  understand 
the  build-up  or  propagation  of  tolerances  and  optimize  the  layout,  features,  or  assembly 
realizations.  This  will  ensure  ease  of  tolerance  delivery.  Any  pro-active  approach  to  assembly 
or  tolerance  analysis  in  the  early  design  stages  will  involve  decision  making  with  incomplete 
information  models.  In  order  to  carry  out  early  tolerance  synthesis  and  assembly  analysis 
in  the  conceptual  design  stage,  we  need  to  devise  a means  of  representation  of  FAB  data 
models  that  will  allow  analysis,  even  with  the  incomplete  data  set. 

The  report  is  organized  as  follows.  In  order  to  build  the  subject  matter  and  explain 
the  FAB  and  Design  For  Tolerancing  (DFT)  models  we  first  explain  the  basic  concepts 
and  methods  of  tolerancing  analysis  and  synthesis  and  describe  the  current  and  evolving 
standards  for  statistical  tolerancing  in  Section  2.  The  proposed  information  model  for  FAB 
is  discussed  in  Section  3.  In  Section  4,  we  present  the  integration  of  FAB  and  DFT  models.  In 
Section  5 we  discuss  about  the  current  and  evolving  standards  and  the  contribution  from  the 
proposed  models  (FAB  and  DFT).  Suggestions  for  further  research  is  presented  in  Section  6. 
Finally  we  present  some  conclusions  in  Section  7. 

A brief  tutorial  on  process  capability  indices  is  given  in  Appendix  A.  The  defect  rate  of 
3.4  parts  per  million  (ppm)  and  the  employed  assumptions  and  convention  is  explained  in 
Appendix  B. 


2 Methods  of  Tolerance  Analysis  and  Synthesis 

Tolerance  design  is  a two-stage  process  consisting  of  tolerance  synthesis  and  analysis.  Toler- 
ance synthesis  is  the  process  of  allocating  tolerances  among  the  geometric  parameters  (dimen- 

1We  have  chosen  the  name  FAB  model  to  draw  the  readers  attention  to  the  fact  that  the  product  design 
activity  starts  with  the  functional  requirements,  assembly  (either  a skeletal  or  list  of  parts),  followed  by  a 
quick  check  on  meeting  the  functional  requirements  by  doing  a behavioral  simulation 
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sions)  of  the  part/component  such  that  the  cost  incurred  due  to  tolerances  is  minimized  and 
at  the  same  time  the  functional  and  assembly  requirements  are  satisfied.  Tolerance  analysis 
on  the  other  hand,  is  concerned  with  the  aggregate  behavior  of  the  individual  tolerances. 

2.1  Methods  for  Tolerance  Analysis 

Tolerance  analysis  can  be  either  worst- case  or  statistical.  Worst-case  tolerance  analysis-also 
called  deterministic  or  high-low  tolerance  analysis-considers  the  worst  possible  combinations 
of  individual  tolerances  and  examines  the  assemblability  of  the  parts,  so  as  to  achieve  100% 
interchangeability  of  parts  in  an  assembly.  This  results  in  unnecessarily  tight  part  tolerances 
and  hence  high  production  costs.  Statistical  tolerancing  is  a more  practical  and  economical 
way  of  looking  at  tolerances  and  works  on  setting  the  tolerances  so  as  to  assure  a desired 
yield.  Here,  the  designer  abandons  the  notion  of  100%  interchangeability  and  accepts  some 
small  percent  of  non-conformance. 

In  tolerance  analysis,  we  seek  to  ascertain  the  validity  of  the  specified  part  tolerances 
with  regard  to  assemblability  and  manufacturability.  Tolerance  analysis  can  be  carried  out 
by  determining  the  tolerance  zones  [65,66]  belonging  to  the  specified  tolerance  types.  It 
can  be  either  statistical  tolerance  analysis  (where  statistical  methods  are  used  together  with 
accompanying  probability  distribution)  or  worst  case  analysis  (i.e.  the  study  of  extreme 
cases).  The  FAB  data  model  as  discussed  in  Section  3.3.1  (page:  38)  provides  all  the  required 
information  for  tolerance  analysis  as  the  data  model  can  be  viewed  as  both  a functional  and 
a structural  hierarchy.  In  order  to  carry  out  tolerance  analysis  in  the  early  design  stages 
when  geometry  is  ill  defined  (i.e.,  incomplete)  the  functional  hierarchy  becomes  very  helpful. 

The  tolerance  analysis  of  the  toleranced  parts  in  assembly  shows  how  a structure  of 
the  part-assembly  satisfies  the  constraints  of  design  and  manufacture,  and  how  the  position 
of  parts  or  their  features  are  constrained  by  the  tolerance  specification.  Tolerance  models 
based  on  well-defined  variational  models  are  needed  to  help  analyze  the  assembly  configu- 
ration uncertainties.  Mathematical  formalization  of  assembly  conditions  is  needed  to  show 
how  contacts  between  tolerance  parts  are  modeled.  The  FAB  data  model  provides  explicit 
descriptions  of  functions  of  the  various  components  and  their  spatial  geometric  relationships, 
including  the  tolerancing  information. 

Statistical  tolerance  analysis  uses  a relationship  of  the  form: 

Y = f{X1 Xn) 

where  Y is  the  response  (a  measurable  characteristic  such  as  assembly  gap)  of  the  assembly 
and  Xi, , Xn  are  the  values  of  some  characteristics  (such  as  dimensions)  of  the  individual 
parts  or  sub-assemblies  making  up  the  assembly.  We  call  / the  assembly  response  function 
(ARF).  The  relationship  can  exist  in  any  form  for  which  it  is  possible  to  compute  a value 
for  Y given  values  of  Xi, . . . , Xn.  It  could  be  any  of  the  following:  an  explicit  or  an  implicit 
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analytic  expression;  a complex  engineering  calculations;  result  from  experiments  and/or 
simulations.  The  input  variables  Xi, . . . ,Xn  are  continuous  random  variables.  In  general, 
they  could  be  mutually  dependent.  The  function  / is  a deterministic  relationship;  Y is  a 
continuous  random  variable.  The  general  problem  of  tolerance  analysis  is  to  compute  the 
probability  distribution  of  Y given  the  distributions  of  Xi, . . . , Xn.  However,  more  often 
we  are  usually  interested  in  computing  the  first  few  moments  (mean,  standard  deviation, 
skewness,  and  kurtosis),  given  the  distributions  or  first  few  moments  of  the  input  variables. 
Once  the  moments  of  Y are  determined,  one  can  compute  a tolerance  range  for  Y that  would 
envelope  a given  fraction  of  the  assembly  yield. 

There  are  a variety  of  methods  and  techniques  available  for  the  above  computational 
problem.  Essentially,  the  methods  can  be  categorized  into  four  classes  [18]: 

1.  Stack  Tolerancing  or  Linear  Propagation  (Root  Sum  of  Squares) 

2.  Non-linear  propagation  (Extended  Taylor  series) 

3.  Numerical  integration  (Quadrature  technique) 

4.  Monte  Carlo  simulation 


2.1.1  Linear  Propagation 

This  is  also  called  stack  tolerancing  and  uses  the  well-known  root  sum  of  squares  (RSS) 
formula.  The  assembly  response  function  here  is  of  the  form: 

Y = a0  + d\Xi  + Q2X2  + . . . + dnXn 

where  a0,  ai, . . . , an  are  constants  and  Xi, . . . , Xn  are  assumed  to  be  mutually  independent. 
Many  dimensional  and  gap-related  measures  fall  into  this  category.  Because  of  the  linear 
relationship  and  mutual  independence,  the  mean  and  variance  of  Y are  given  by: 


/iy  — do  + d\ /J,i  + d2^2  + • ■ • + dn[ln 

2 2 2 1 2 2 1 ,22 

(Jy  = al(Tl  + a2cr2  + . . . + dnan 


where  / ii  is  the  mean  and  cq,  the  standard  deviation  of  Xi,  i = 1 , ,n.  The  nomenclature 
RSS  arises  because  of  the  above  formula  for  standard  deviation.  If  the  individual  distribu- 
tions are  normal,  then  Y is  also  normally  distributed.  Even  if  the  individual  distributions 
are  not  normal,  Y can  safely  be  treated  as  normal,  by  invoking  the  central  limit  theorem. 

If  the  linear  relation  for  Y above  is  only  approximately  true,  then  one  can  expand 
f(Xi, . . . , Xn)  as  a Taylor  series  and  drop  all  but  the  constant  and  linear  terms.  This  is 
often-used  device  in  statistical  tolerancing  to  handle  approximately  linear  relationships.  In 
such  a case, 


evaluated  at  aq 


liii  i — 1, . . . , n, 
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and  all  of  the  constant  terms  are  gathered  into  ao-  The  computation  of  the  above  partial 
derivatives  could  be  of  two  types.  In  the  first  case,  the  function  / is  known  and  the  partial 
derivatives  are  known  to  exist.  In  the  second  case,  the  functional  relationship  is  either  too 
intractable  or  not  even  available  in  analytic  form.  In  such  a case,  numerical  estimates  have 
to  be  obtained  for  the  partial  derivatives  [18]. 

The  linear  case  is  the  simplest  and  the  most  efficient  among  all  tolerance  analysis  ap- 
proaches. It  is  very  appealing  for  synthesis  methods  that  use  analysis  in  an  iterative  way. 


2.1.2  Non-linear  Propagation  (Extended  Taylor  Series) 

If  the  assembly  response  function  Y is  highly  non-linear,  application  of  the  RSS  method 
could  lead  to  serious  errors.  An  extended  Taylor  series  approximation  for  the  relationship 
/ can  possibly  be  employed  provided  / is  available  in  an  analytic  form.  Usually,  the  ex- 
pansion is  considered  up  to  the  sixth  order.  The  expansion  is  possible  only  when  all  the 
appropriate  partial  derivatives  exist.  The  main  computational  issue  involves  computing  the 
partial  derivatives  of  /.  Mathematically  tractable  formulae  for  the  first  four  moments  of  Y 
are  available  [18]  and  are  ideally  suited  for  tolerance  analysis  and  synthesis.  These  formulae 
need  only  the  first  four  moments  of  the  distributions  of  Xi, ...  ,Xn.  Most  often,  the  partial 
derivatives  are  computed  using  analytic  methods.  If  closed  form  solutions  do  not  exist  then 
we  need  to  use  numerical  iterations  techniques  as  described  in  the  next  section. 


2.1.3  Numerical  Integration 

If  the  function  / is  not  available  in  analytic  form  and  Y can  only  be  computed  through 
numerical  calculations  or  engineering  methods  or  simulations,  numerical  methods  have  to  be 
used.  Quadrature  methods  are  generally  used  here.  Here,  we  assume  that  for  any  function 
h(Xi, . . . , Xn)  (different  from  /)  of  mutually  independent  random  variables  Xi, . . . , Xn  with 
probability  density  functions  wxt(xi),  the  expected  value  of  h is  given  by  the  integral 


71 

) E[  (' wxi(xi)dxi ) 

i= 1 


The  above  expression  can  be  approximated  by  a quadrature  expression  [18]  that  involves 
evaluations  of  h at  2 n2  + 1 prescribed  values.  These  evaluations  involve  only  the  first  four 
moments  of  X\, . . . ,Xn.  Given  an  assembly  response  function  /,  a corresponding  function 
h as  above  can  be  defined  and  simple  moment  transfer  relations  can  be  used  to  compute 
the  first  four  moments  of  /.  The  quadrature  technique  adapts  well  to  statistical  tolerancing 
problems  since  it  can  handle  the  iteration  inherent  in  a tolerancing  problem  efficiently.  An 
improved  integration  technique  is  provided  in  [55,58]. 
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2.1.4  Monte  Carlo  Simulation 

The  appeal  of  Monte  Carlo  simulation  lies  in  its  applicability  under  very  general  settings 
and  the  unlimited  precision  that  can  be  achieved.  In  particular,  Monte  Carlo  simulation  can 
be  used  in  all  situations  in  which  the  above  three  techniques  (stack  tolerancing,  extended 
Taylor  series,  and  numerical  integration)  can  be  used  and  can  yield  more  precise  estimates. 
For  this  reason,  the  Monte  Carlo  technique  is  easily  the  most  popular  tool  used  in  toleranc- 
ing problems.  The  caveat,  however,  is  the  large  computational  time:  for  situations  where 
the  above  three  techniques  are  adequate  and  have  acceptable  precision,  the  Monte  Carlo 
technique  is  much  more  expensive  in  terms  of  computational  time. 

Monte  Carlo  simulation  proceeds  as  follows  (Figure  1).  Pseudo  random  number  gen- 
erators are  used  to  generate  a sample  set  of  numbers  Xi, ...  ,xn,  belonging  to  the  random 
variables  Xi, . . . ,Xn,  respectively.  The  value  of  Y,  say  yi  = /(aq, . . . ,xn ),  corresponding 
to  this  sample  is  computed.  This  procedure  is  replicated  a large  number  of  times,  say  N 
times.  This  would  yield  a random  sample,  {yi, . . . ,yv}i  for  Y.  Standard  statistical  estima- 
tion methods  are  then  used  to  analyze  the  distribution  of  Y.  The  precision  of  this  statistical 
analysis  increases  as  proportional  to  VN  and  therefore  unlimited  precision  can  be  achieved 
through  large  number  of  replications.  Special  techniques  are  available  for  significantly  en- 
hancing the  precision  of  the  Monte  Carlo  method  for  a given  N.  These  include:  weighted 
sampling,  reuse  of  samples,  and  use  of  approximation  functions  [18]. 


Figure  1:  Flow  chart  of  Monte-Carlo  simulation  method 
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2.2  Methods  for  Tolerance  Synthesis 

In  the  context  of  electro-mechanical  assembly  design,  tolerance  synthesis  usually  refers  to  the 
allocation  of  specified  assembly  tolerances  among  the  constituent  parts  and  sub-assemblies, 
so  as  to  ensure  a specified  yield  or  minimize  a proper  cost  function.  More  generally,  if 
Y = is  an  assembly  response  function,  then  the  synthesis  problem  involves 

finding  the  best  nominals  and  tolerances  for  Xi, . . . , Xn , given  nominal  and  tolerance  specifi- 
cations for  Y . Synthesis  is  naturally  an  optimization  problem;  one  can  formulate  an  objective 
function  that  captures  yield  requirements  or  production  cost  requirements  and  pose  an  op- 
timization problem  by  including  tolerance  related  constraints. 

There  are  several  views  and  variants  of  the  synthesis  problem,  depending  on  the  objective 
function  and  the  constraints.  One  view  is  to  minimize  the  total  manufacturing  cost  by 
choosing  both  the  individual  part  tolerances  and  the  manufacturing  processes  for  making 
the  individual  parts.  This  requires  cost  versus  tolerance  relationships  for  each  individual 
dimension.  Another  view7  is  to  find  robust  nominals  for  individual  dimensions,  i.e. , nominal 
values  at  which  the  effect  of  variations  on  the  assembly  response  function  is  minimum.  This 
issue  is  addressed  by  Taguchi’s  robust  design  methodology  and  Park's  Holistic  Probabilistic 
Design  (HPD).  Also,  depending  on  the  nature  of  the  objective  function  and  the  constraints, 
the  synthesis  problem  can  be  deterministic  or  stochastic. 

To  formulate  the  synthesis  problem  meaningfully,  a certain  amount  of  preprocessing 
is  often  required.  For  example,  one  needs  to  first  determine  the  tolerance  limits  on  the 
assembly  response  function,  Y.  An  important  preprocessing  step  is  sensitivity  analysis , which 
determines  which  assembly  parameter  variations  have  significant  effects  on  the  assembly 
response  function.  This  reveals  the  set  of  parameters  or  individual  dimensions  to  emphasize 
in  the  synthesis  procedure. 

2.2.1  Iterative  Methods  Based  on  Analysis 

A simple  and  realistic  mechanism  for  tolerance  synthesis  is  to  employ  a trial  and  error  tech- 
nique for  postulating  tolerances  for  individual  parts  and  sub-assemblies.  Next,  perform  a 
statistical  tolerance  analysis  required  to  ascertain  whether  this  postulated  set  of  tolerances 
fulfills  the  desired  criteria.  If  the  chosen  set  is  unsatisfactory,  some  of  the  tolerances  are 
changed  and  the  analysis  redone;  this  step  is  repeated  until  a satisfactory  set  of  part  toler- 
ances is  obtained.  Typically,  at  the  end  of  each  iteration,  we  obtain  a probability  of  assembly, 
a probability  of  conformance,  an  expected  yield,  or  a more  detailed  cost.  This  technique  ef- 
fects trade-offs  during  each  iteration  and  is  appealing  as  it  uses  the  findings  of  the  current 
iteration  to  drive  the  next  iteration.  During  early  iterations,  approximate  cost  figures  and 
less  accurate  estimates  can  be  used.  These  can  be  replaced  by  more  accurate  figures  as  the 
iterations  start  producing  good  solutions. 

The  methods  discussed  for  statistical  tolerance  analysis,  namely  stack  tolerancing,  ex- 
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tended  Taylor  series,  quadrature  methods,  and  Monte  Carlo,  are  all  suited  for  the  iterative 
approach.  Evans  [18]  has  discussed  the  merits  and  issues  concerning  the  use  of  these  methods 
for  the  iterative  methodology  for  synthesis. 

2.2.2  Multistage  Tolerance  Synthesis 

In  order  to  synthesize  tolerance  in  the  conceptual  design  stage,  Roy  and  Bharadwaj  [63] 
suggested  a multi  stage  procedure.  Figure  2 shows  the  conceptual  schema  for  tolerance 
synthesis.  Given  the  design  function  requirements,  manufacturing  processing  information 
and  assembly  plan,  the  schema  helps  assign  both  dimensional  and  geometric  tolerances  (along 
with  required  datum  reference  planes)  to  be  part  of  an  assembly.  The  tolerance  synthesis 
schema  starts  with  collecting  the  following  information  (refer  to  Phase  1,  Figure  2)  from  the 
FAB  data  model. 

1.  Geometry  description:  the  schema  requires  geometry  description  at  the  following  levels: 

• assembly  - position  and  orientation  information  for  each  component  within  the 
assembly 

• part-spatial  location  of  form  features  in  the  part  and  their  inter-relationships  and 

• feature-feature  geometry. 

2.  Part  function  specifications:  Functional  and  behavioral  attributes  of  each  part  in  the 
assembly. 

3.  Material  and  surface  finish  specifications: 

• Material  and  surface  characteristics  should  be  either  retrieved  from  the  database 
or  supplied  manually  by  the  user. 

4.  Assembly  graph: 

• The  procedure  for  assembling  different  parts  in  the  assembly  (without  considering 
the  effect  of  tolerances)  should  be  retrieved  from  the  assembly  data  model.  The 
assembly  graph  denotes  the  assembly  sequence  of  the  parts  in  order  to  form  an 
assembly.  It  is  a directed  acyclic  graph  whose  nodes  are  represented  by  “assembly 
states”  and  whose  edges  denote  the  ’’assembly  process”. 

In  the  second  phase,  the  part  function  model  (PFM)  and  the  process  model  for  each 
part  of  the  assembly  are  generated.  A part  function  model  describes  the  spatial  and  design 
relationships  that  exist  on  the  mating  faces  of  a part  in  terms  of  certain  kinematics,  force 
degrees  of  freedom  (dofs),  presence/absence  of  motions,  and  transmission  of  forces  along 
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PHASE  1 SPECIFICATION 


PHASE  2 GENERATION 


PHASE  3 SYNTHESIS 


PHASE  4 ASSIGNMENT 


Figure  2:  Tolerance  Synthesis  Scheme 


the  particular  axes  of  a surface  [7].  The  process  model  represents  the  process  plan  for 
manufacturing  the  part  without  considering  the  effect  of  tolerance  [4], 

The  third  phase  of  the  schema  is  the  synthesis  stage.  Different  types  of  tolerances  are 
synthesized  for  each  part  of  the  assembly.  It  consists  of  two  major  tasks:  (i)  transforma- 
tion of  the  PFM  model  into  functional  tolerance  limits,  and  (ii)  constraining  the  functional 
tolerance  limits  with  respect  to  different  manufacturability  and  assemblability  constraints. 
The  first  task  can  be  achieved  by  developing  appropriate  application  domain-specific  PFM- 
to-Functional-Tolerance-Limit  maps  (refer  to  [63,64]  for  a detailed  discussion).  The  second 
task  can  be  achieved  by  developing  optimization  problems  which  contain  both  the  functional 
tolerance  limits  and  the  different  constraints.  In  the  fourth  phase,  dimensional  and  geomet- 
ric tolerances  (along  with  the  datum  specifications)  are  finetuned  with  respect  to  the  design 
functions  and  manufacturing  constraints. 

The  FAB  data  model  provides  the  required  information  needed  for  the  above  four  stage 
tolerance  synthesis  process. 
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2.2.3  Optimization  Methods 


Since  tolerance  synthesis  can  be  posed  as  an  optimization  problem,  mathematical  program- 
ming techniques  such  as  linear  programming,  non-linear  programming,  and,  integer  pro- 
gramming are  relevant.  There  have  been  several  efforts  in  this  direction  [13,43,46].  Also, 
heuristic  techniques  for  optimization  such  as  simulated  annealing,  genetic  algorithms,  La- 
grangian  relaxation,  and  Tabu  search  have  been  used  by  researchers  [25,38,45]. 

The  following  integer  programming  formulation  given  by  Kusiak  and  Feng  [43],  provides 
a flavor  of  a typical  synthesis  problem  in  the  optimization  framework: 

Minimize 


subject  to: 


E E CijXij 
i=lj=l 


y y Uj%ij  < t 

i=l  j= 1 


E xij  = 1 

3 = 1 

Xij  = 0,1  \/l:j 

where  the  index  i denotes  one  of  n dimensions  involved  in  the  assembly;  index  j denotes 
one  of  m manufacturing  processes  that  can  be  used  for  producing  a dimension;  Cy  is  the 
manufacturing  cost  of  process  j when  used  for  producing  dimension  z;  Uj  is  the  3a  normal 
variation  of  process  j when  used  to  produce  dimension  z;  T is  the  tolerance  stack-up  limit 
for  the  assembly;  and  xi3  is  a binary  decision  variable  that  takes  a value  1 if  process  j is 
selected  for  producing  dimension  z and  0 otherwise.  Note  that  the  objective  is  to  minimize 
the  total  direct  manufacturing  cost,  by  choosing  the  appropriate  tolerances  and  the  right 
mix  of  manufacturing  processes.  In  the  above  formulation,  a linear  relationship  has  been 
assumed  between  part  tolerances  and  also  worst  case  tolerancing  has  been  used.  Hence,  it  is 
a deterministic  tolerance  synthesis  problem.  In  HPD,  the  optimization  is  done  considering 
nominals  and  variance  together. 


2.2.4  Design  of  Experiments 

Here,  the  assembly  response  function  (or  in  general,  a well-defined  cost  function),  is  com- 
puted for  various  discretized  values  of  the  random  variables  Xi, . . . , Xn,  (dimensions  with 
tolerances)  according  to  the  design  of  experiments  theory.  The  factors  used  in  the  experiment 
include  not  only  the  individual  values  of  Xi, . , . , Xn , but  also  factors  that  capture  tolerance 
related  constraints.  Full  factorial  or  fractional  factorial  designs  can  be  used  depending  on 
the  number  of  factors  and  levels  of  the  factors.  Prior  sensitivity  analysis  can  help  in  choosing 
the  levels  for  the  factors.  The  setting  that  leads  to  a minimum  cost  and  also  satisfies  the 
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tolerance  constraints  can  be  chosen  as  the  solution.  One  can  go  a step  further  and  fit  a 
statistical  model  that  describes  the  cost  function  in  terms  of  all  the  factors  and  use  this 
model  to  arrive  at  an  optimal  solution  for  the  problem. 

Taguchi  methods,  which  are  described  in  the  next  section,  use  design  of  experiments  in 
a novel  way  to  find  robust  nominals. 

2.3  Best  Practices 

In  the  last  decade,  many  companies  have  established  comprehensive  programs  in  total  quality 
management.  These  efforts  include  those  of  Motorola,  Xerox,  IBM,  AT&T  Bell  Laborato- 
ries, and  several  others  that  have  initiated  formal,  corporate  programs  for  improved  tolerance 
specification,  monitoring,  and  control.  In  this  section,  we  outline  the  tolerancing  best  prac- 
tices developed  at  Motorola  and  Xerox  which  are  closely  related  to  our  approach.  We  also 
provide  an  overview  of  Taguchi’s  robust  design  methodology,  which  has  emerged  as  a best 
practice  in  many  companies. 

2.3.1  Motorola  Six  Sigma  Program 

Six  sigma  quality  is  the  benchmark  of  excellence  for  product  and  process  quality,  popularized 
by  Motorola  [22,24].  It  provides  a quantitative  and  statistical  notion  of  quality  useful  in 
understanding,  measuring,  and  reducing  variation.  A product  is  said  to  be  of  six  sigma 
quality  if  there  are  no  more  than  3.4  non-conformities  per  million  opportunities  (3.4  ppm) 
at  the  part  and  process-step  level,  in  the  presence  of  typical  sources  of  variation.  The  six 
sigma  quality  concept  recognizes  that  variations  are  inevitable  due  to  insufficient  design 
margin,  inadequate  process  control,  imperfect  parts,  imperfect  materials,  fluctuations  in 
environmental  conditions,  operator  variations,  etc. 

Tolerance  analysis  and  synthesis  in  the  six  sigma  program  are  based  on  the  following 
criteria: 

1.  six  sigma  characterization  of  products  and  processes  where  the  process  capability  in- 
dices Cp  and  CPk  are  used  as  the  vehicles  to  characterize  the  product-process  quality; 

2.  simple,  intuitive  extensions  to  the  RSS  method  which  enables  tolerance  analysis  and 
synthesis  in  the  presence  of  shifts  and  drifts  of  the  process  mean;  and 

3.  a well-defined,  systematic  program  for  design  for  quality  taking  into  account  both  the 
product  perspective  and  the  process  perspective. 

We  provide  a brief  outline  of  these  issues;  for  a detailed  discussion  refer  [22-24].  Before 
discussing  the  tolerance  analysis  and  synthesis  using  Motorola  six-sigma  methodology,  we 
discuss  the  Motorola  six-sigma  quality  program. 
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Definition  of  Motorola  Six-Sigma  Quality 

According  to  the  Motorola  six-sigma  program  [22],  a product  is  said  to  be  of  six-sigma 
quality  (0.9999966  confidence)  if  there  are  no  more  than  3.4  non-conformities  per  million  of 
opportunities,  at  the  product  and  process  level,  in  the  presence  of  typical  sources  of  variation 
(shifts  and  drifts).  The  six-sigma  quality  represents  a 2941  (10000/3.4)  fold  improvement  in 
quality  compared  to  the  99%  quality. 

The  six-sigma  design  uses  the  capability  indices  Cp  and  Cpk  as  quality  metrics  (refer 
Appendix  A,  page  67  for  more  details).  The  six-sigma  design  has  a defect  rate  of  3.4  parts 
per  million  with  Cp  > 2 and  CPk  >1.5  under  the  following  assumptions  that  the  process 
(or  a typical  product  characteristic)  (1)  is  normally  distributed  and  (2)  has  a mean  shift  of 
1.5(7  units  from  the  target  value.  The  explanation  of  the  defect  rate  3.4  ppm  with  aid  of  an 
example  is  given  in  Appendix  B. 

2.3.2  Tolerance  Analysis  and  Synthesis  in  the  Motorola  Approach 

The  Motorola  program  assumes  a linear  model  for  F,  of  the  form 


1'  — <2o  + aiJAi  + 0,2X2  T- . • . + onXn 

If  there  is  no  mean  shift,  then  the  standard  RSS  formulae  are  applicable: 


/iy  — Oo  + ai/Ii  + a2/i2  + • • • + On/Oji 

a2Y  = a\o\  + a\a\  + . . . + 

Recall  that  cq,  for  i — 1 , ...  ,n,  can  also  be  written  as: 


where  Tj  is  the  tolerance  range  of  the  ith  part  and  Cpi  is  the  Cp  value  for  the  zth  part 
(i  = 1, ...  ,n).  In  the  presence  of  a mean  shift,  the  standard  RSS  formula  cannot  be  used 
for  computing  the  standard  deviation  of  Y.  Two  alternative  approaches  are  recommended 
by  the  Motorola  program.  The  first  is  the  Dynamic  RSS  (DRSS)  where  the  CPk  values, 
CPki , . . . , CPkn , of  the  individual  processes  corresponding  to  dimensions  Xi,...,Xn  and  the 
tolerances,  T-L, ...  ,Tn,  of  the  individual  parts  are  used  in  the  following  way  to  compute  the 
variance  of  Y: 
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Note  that  the  standard  deviations  cq  are  amplified  by  an  amount  equal  to  , for  i = 
1, . . . , n.  Thus  the  DRSS  method  emulates  random  behavior  in  the  process  mean  by  inflating 
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the  process  standard  deviation,  but  has  little  impact  on  the  overall  mean.  The  second 
alternative  method  called  Static  RSS  (SRSS)  achieves  the  sustained  mean-shifts  condition 
by  applying  a correction  factor  to  the  individual  nominals.  For  details,  see  [24]. 

Tolerance  analysis  is  carried  out  by  using  RSS,  DRSS,  and  SRSS,  as  appropriate.  Tol- 
erance synthesis  uses  tolerance  analysis  in  an  iterative  way.  Each  iteration  will  evaluate 
the  resulting  probability  of  non-conformance  and  the  Cp  and  Cpk  values.  The  goal  of  the 
synthesis  procedure  is  to  obtain  a probability  of  non-conformance  of  at  most  3.4  ppm,  which 
is  guaranteed  by  Cp  > 2 and  Cpk  >1.5.  The  synthesis  can  assume  several  forms: 

1.  finding  optimal  values  for  nominal  dimensions; 

2.  finding  optimal  values  for  tolerances;  and 

3.  establishing  a variance  pool  that  can  be  allocated  to  individual  processes  so  as  to  obtain 
the  desired  assembly  yield. 

The  Motorola  six  sigma  approach  uses  the  normal  distribution  for  all  its  probability  and 
tolerancing  computations.  While  this  can  be  listed  as  a limitation,  it  takes  very  little  away 
from  the  intrinsic  novelty  and  applicability  of  the  approach.  The  ideas  it  has  germinated 
essentially  hold  in  all  situations;  only  the  probability  computations  need  to  be  redone  under 
non- normal  situations  and  the  quantitative  measures  need  appropriate  reinterpretation. 

2.3.3  Holistic  Probabilistic  Design  (HPD) 

The  HPD  methodology  [56-58]  is  one  of  several  quality  programs  at  the  Xerox  Corporation. 
The  program  is  based  on  relating  service  dissatisfiers  and  customer  tolerances  to  a set  of 
critical  parameters  (parameters  that  are  critical  to  the  product’s  function).  The  tolerances 
of  the  critical  parameters  are  related  to  piece  part  variabilities  through  multiple  flow  of 
variance  chains.  Tolerance  analysis  and  synthesis  are  carried  out  through  the  chains  to  yield 
the  desired  quality.  Since  the  objective  is  to  maximize  the  amount  of  manufacturing  and 
usage  variability  the  product  can  tolerate,  with  negligible  impact  on  the  targeted  level  of 
performance,  the  program  is  also  called  design  for  latitude.  The  methodology  is  implemented 
using  a complete  suite  of  tools  for  stochastic  variability  analysis.  These  include  the  following: 

1.  A stochastic  modeling  based  technique  for  computing  the  distribution  of  a function  of 
a random  variables  almost  exactly; 

2.  Contribution  analysis  that  provides  a reliable  guidance  on  factors  that  have  a significant 
effect  on  the  assembly  response  function;  and 

3.  An  operating  window  optimization  method  that  helps  choose  the  densities  of  certain 
input  random  variables  for  which  the  allowable  range  of  operation  is  maximized. 
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Tolerance  analysis  is  based  on  a stochastic  technique  that  uses  a failure  rate  prediction 
model.  Let  Y — f(Xi,...,  Xn)  be  the  assembly  response  function,  as  usual.  If  failure  is 
defined  by  /(Xl, . . . , Xn)  > YQ,  then  the  probability  of  failure  is  given  by: 

Pr{f{X1,...,Xn)>Y0}  = /.../  w{X1,...,Xn)dX1...dXn 

J a.  J 

where  D is  the  n-dimensional  failure  region  and  w(X i, . . . , Xn)  is  the  joint  density  of  the  n 
random  variables.  From  this,  it  is  easy  to  see  how  to  compute  the  distribution  function  of 
Y.  Assuming  mutual  independence  of  Xi, . . . , Xn,  the  HPD  tool  uses  an  efficient  numerical 
technique  to  evaluate  such  multiple  integrals.  This  enables  us  to  compute  the  distribution 
accurately.  This  computation  is  versatile  since  it  can  handle  any  type  of  distribution  and  any 
type  of  relationship;  it  has  excellent  computational  performance  if  the  number  of  random 
variables  is  less  than  10.  The  above  computation  enables  variability  analysis,  hence  tolerance 
analysis.  It  also  provides  a sound  basis  for  iterative  tolerance  synthesis.  An  attractive  feature 
of  this  technique  is  its  applicability  to  both  geometric  and  non-geometric  type  of  situations. 
For  example,  the  random  variables  Xi, ...  ,Xn  need  not  be  dimensions;  they  could  be  physical 
quantities  such  as  force,  pressure,  temperature,  and  speed. 

The  contribution  analysis  embedded  within  HPD  is  a powerful  feature  of  HPD.  It  pro- 
vides a sound  basis  for  determining  the  input  variables  that  have  a pronounced  effect  on 
the  assembly  response  function.  Also,  it  takes  into  account  the  nature  of  the  input  distri- 
butions and  accounts  for  cross-term  effects.  This  has  several  advantages  over  other  existing 
techniques  for  contribution  analysis. 

Another  feature  of  HPD  is  the  stochastic  operating  window  optimization.  This  feature 
enables  us  to  maximize  the  allowable  range  of  operation  by  intelligently  selecting  the  densities 
of  input  random  variables  using  tools  provided  by  HPD. 

The  HPD  tool  consists  of  two  major  modules:  HPD-VA  and  HPD-OPT.  The  module 
HPD-VA  is  a stochastic  analyzer  that  includes  variability  analysis  and  tolerance  analysis. 
HPD-OPT  has  a wide  variety  of  deterministic,  stochastic,  and  statistical  optimization  rou- 
tines. HPD-OPT  finds  the  most  robust  set  of  nominals  and  tolerances. 

2.3.4  Taguchi  Methods 

Taguchi  methods,  also  known  as  robust  design  methods,  [36,37,59],  are  technical  methods 
for  quality  and  cost  control  at  the  product  and  process  design  stages.  According  to  Taguchi, 
the  cost  of  a product  is  the  loss  incurred  to  the  society  before  the  product  is  shipped  to 
the  customer,  whereas  quality  is  the  loss  imparted  to  the  society  after  the  product  has  been 
shipped  to  the  customer.  Such  losses  include  the  following:  loss  due  to  harmful  side  effects; 
loss  due  to  variations  in  the  product’s  performance  characteristics;  and  all  losses  that  can 
be  traced  to  the  poor  performance  of  the  product.  Taguchi  methods  emphasize  reducing 
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the  sensitivity  of  engineering  designs  to  various  sources  of  variation.  The  methods  are  cost- 
effective  as  it  minimizes  the  influence  of  variation  sources  rather  than  control  them. 

Let  Y be  a performance  characteristic ; as  before,  Y is  a continuous  random  variable. 
Taguchi  considers  Y as  a function  of  design  parameters  or  design  factors , 0 = (6\, . . . , 0*), 
and  noise  parameters  or  noise  factors , W = (w1, . . . , wt).  Thus,  Y = f{&,W).  Design 
factors  are  input  variables  whose  nominal  settings  and  that  have  a pronounced  influence  on 
Y can  be  chosen  by  a designer.  Design  factors  are  of  two  types:  signal  factors , which  affect 
only  the  mean  of  Y,  and  control  factors , which  affect  both  mean  and  variance  of  Y.  Noise 
factors  are  input  variables  that  cause  Y to  deviate  from  its  target  value.  Noise  factors  include 
deviations  of  the  actual  values  of  design  factors  from  the  nominal  settings.  Taguchi  considers 
two  types  of  matrices:  the  design  matrix  (inner  array)  and  the  noise  matrix  (outer  array). 
The  design  matrix  has  k columns,  each  column  corresponding  to  a particular  design  factor. 
Each  row  of  this  matrix  represents  a specific  combination  of  design  parameter  settings.  The 
number  of  rows  depends  on  the  number  of  combinations  of  design  parameter  settings  being 
investigated.  Similarly,  the  noise  matrix  has  t columns,  each  column  corresponding  to  a 
particular  noise  factor.  Each  row  of  this  matrix  represents  a specific  combination  of  noise 
factor  settings.  The  number  of  rows  depends  on  the  number  of  combinations  of  noise  factor 
settings  sought  to  be  investigated. 

Let  r be  the  target  value  (ideal  value)  of  Y;  p,  its  mean;  and  cr,  its  standard  deviation. 
The  target  value  need  not  be  the  midpoint  of  a tolerance  interval.  Variations  of  Y about  the 
target  value  r cause  losses  to  the  customers.  Let  l(Y)  be  the  loss  due  to  the  deviation  of  Y 
from  r.  Taguchi  suggests  a quadratic  form  for  the  loss  function: 

l(Y)  = k(Y  — t)2 

where  A:  is  a constant  that  can  be  computed  from  a known  value  of  loss  at  any  designated 
value  of  Y.  The  loss  function  is  a random  variable  and  the  expected  quadratic  loss,  E[l(Y)], 
is  given  by 

L = E[l(Y)}  = k(o2  + {p-  r)2) 

Thus  the  expected  quadratic  loss  is  the  the  sum  of  the  variance  of  Y and  the  square  of  bias 
(bias  indicates  how  far  away  from  the  target  value  the  process  mean  is  behaving). 

Minimization  of  the  expected  quadratic  loss  is  the  primary  objective  of  Taguchi  methods. 
This  is  accomplished  by  maximizing  a signal-to-noise  ratio  (also  called  as  a performance 
statistic).  Taguchi’s  use  of  signal-to-noise  ratios  rather  than  directly  employing  the  expected 
quadratic  loss  is  motivated  by  considerations  such  as  ease  of  statistical  estimation,  more 
direct  coupling  to  design  factors,  and  improved  additivity  of  control  factor  effects.  For  a 
detailed  exposition  refer  [37,  59].  A signal-to-noise  ratio  is  a statistical  estimate  of  the  effect 
of  noise  factors  on  Y for  a particular  setting  of  design  parameters.  Numerous  performance 
statistics  have  been  defined  by  Taguchi  (more  than  60). 

The  following  are  the  main  steps  in  the  Taguchi  method. 
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1.  Identify  appropriate  loss  function  or  signal-to-noise  ratio,  initial  and  competing  settings 
of  design  factors,  and  important  noise  factors  and  their  ranges. 

2.  Construct  the  design  matrix  and  the  noise  matrix.  The  design  matrix  is  chosen  based 
on  the  theory  of  design  of  experiments  or  is  chosen  from  Taguchi’s  collection  of  or- 
thogonal arrays  [37].  The  noise  matrix  is  usually  chosen  from  Taguchi’s  collection  of 
orthogonal  arrays. 

3.  Conduct  a parameter  design  experiment.  This  involves  Nd  runs , where  Nd  is  the 
number  of  rows  of  the  design  matrix.  Each  run  corresponds  to  a particular  row  and 
involves  Nn  replications,  where  Nn  is  the  number  of  rows  of  the  noise  matrix.  For  each 
run,f,  {i  — 1, . . . , Nd ),  a corresponding  signal-to-noise  ratio,  [Z(0)]j  is  computed. 

4.  Use  [^(0)^, . . . , [ Z(0)]jvd , to  predict  a statistical  model  for  the  signal  to  noise  ratio. 
Use  the  predicted  statistical  model  to  obtain  optimal  or  best  design  parameter  settings: 

0*  = .,0*) 

5.  Conduct  a verification  experiment  to  confirm  that  0*  indeed  minimizes  the  expected 
loss.  Otherwise,  iterate. 

Taguchi  methods  make  several  assumptions.  Examples  of  these  include:  absence  of 
interaction  effects  among  the  factors;  additivity  of  control  factors;  separability  of  signal 
factors  and  control  factors;  and  use  of  signal-to-noise  ratios  instead  of  direct  measures  [37, 
59].  However,  the  methodology  embodies  sound  engineering  considerations  and  intuition  for 
obtaining  robust  designs,  which  explains  its  widespread  use.  From  a tolerancing  viewpoint, 
Taguchi  methods  provide  a valuable  tool  for  synthesizing  robust  nominals.  Also,  the  methods 
can  be  applied  potentially  during  early  assembly  design  stages.  Furthermore,  the  methods 
enable  economic  considerations  to  be  incorporated  into  tolerance  analysis  and  synthesis,  and 
provide  an  approach  that  is  complementary  in  many  ways  to  other  best  practices  such  as 
the  Motorola  six  sigma  program  and  the  Xerox  HPD  methodology. 

2.4  Summary 

In  summary,  in  this  section  we  discussed  the  various  tolerance  analysis  and  synthesis  tech- 
niques, tools  and  industry  best  practices.  In  the  next  section  we  will  discuss  about  the 
current  and  evolving  statistical  tolerancing  standards  syntax  and  semantics. 

Although  statistical  tolerancing  has  been  practiced  in  industry  for  a long  time,  standards 
are  not  yet  available.  ISO  [34,72,73]  is  investigating  how  to  standardize  statistical  toler- 
ancing. The  scope  of  ISO  TC  213  [34]  is  standardization  in  the  field  of  geometrical  product 
specifications  (GPS)  i.e.  macro-  and  micro-geometry  specifications  covering  dimensional  and 
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geometrical  tolerancing,  surface  properties  and  the  related  verification  principles,  measuring 
equipment  and  calibration  requirements  including  the  uncertainty  of  dimensional  and  geo- 
metrical measurements.  This  standardization  includes  the  basic  layout  and  explanation  of 
drawing  indications  (symbols). 

2.5  Statistical  Tolerancing 

Both  classical  and  statistical  tolerancing  are  currently  practiced  in  the  industry.  Classical 
tolerancing  is  popularly  known  as  the  worst-case  tolerancing.  Given  an  actual  feature,  it 
defines  a test  to  decide  whether  the  feature  is  acceptable.  Statistical  tolerancing,  on  the 
other  hand,  deals  with  a population  of  features.  It  defines  a test  that  decides  whether  a 
given  population  of  actual  features  is  acceptable. 

Let  us  consider  the  classic  example  discussed  by  Evans  [17].  In  that  example  a stack  of 
n nominally  identical  disks  is  analyzed.  For  illustrative  purpose,  consider  n — 10  and  the 
height  of  the  stack  as  the  response  (assembly  response  function  Y — X i + X2  + - ■ ■ + Xn  [53]). 
Let  us  assume  that  the  stack  height  must  be  1.25  ± 0.01  ( USLy  — 1.26,  LSL  = 1.24).  The 
worst  case  tolerancing  of  the  disks  would  yield  0.125  ± 0.001  for  the  components  (individual 
disks). 

The  problem  with  statistically  tolerancing  the  components  of  an  assembly  (individual 
disks)  deals  with  setting  limits  on  the  allowed  excursion  of  the  components’  mean  and  range 
(using  the  standard  deviation).  Let  us  assume  that  X{  ~ 7V(/jj,  cq)  with  the  assumption  that 
a i = gc  Vi.  If  erf  is  the  variance  of  the  component  i,  the  variance  of  assembly  is  given  by 

cr  = Var  (assembly)  = ^2  erf  = nerf  if  cq  = erc;  Vi  (1) 

Component  level  to  assembly  Level 

Traditionally,  the  half  width  of  the  specified  tolerance  is  set  to  three  times  the  standard 
deviation  of  the  assembly  (Cp  = 1)  in  determining  the  standard  deviation  of  each  component. 
In  this  example,  let  \ii  — 0.125  and  cq  = pp-  so  that  \xy  — X.  Mi  = 1-25  and  ay  — \/T0  . 

That  is,  given  X{  ~ JV(0.125,  ^j^)  then  the  assembly  variable  Y ~ N(  1.25,  \/l0  (^p))-  We 
can  now  compute  the  following  probability: 

Pr{Y  £ [1.24, 1.26]}  « Pr{1.25  - lOcq,  <Y<  1.25  + 10ay}  « 10"20 

Assembly  level  tolerancing  to  component  level 

Now  let  us  do  the  statistical  analysis  from  the  assembly  level  to  component  level.  Set 
±3a  limits  for  the  stack  height  as  0.01.  Since  erf  = Vrar(assembly)  — erf  — no2,  the 
process  would  be  designed  so  that: 

0.01 
Gc  ~ 3\/fo 
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The  tolerance  on  individual  disks  could  be  3 uc  that  is  in  terms  of  high/low  tolerancing 
0.125  ± 0.003.  Thus  approached  from  the  viewpoint  of  statistical  tolerancing,  the  tolerances 
on  the  individual  disks  can  be  relaxed  by  a factor  of  three. 

There  are  various  reasons  for  doing  statistical  tolerancing:  (1)  worst-case  tolerancing  leads 
to  unnecessary  tight  tolerances,  costing  more  money  and  time;  (2)  rate  of  miniaturization 
exceeds  the  rate  of  improvement  in  the  manufacturing  process  capabilities;  and  (3)  it  is  the 
missing  link  in  Shewhart’s  quality  control  process  [75]. 

2.5.1  Statistical  Tolerancing  - Syntax  and  Semantics 

The  boundary  of  an  actual  part  (model  of  a manufactured  part)  can  be  partitioned  into 
many  actual  surface  features.  Since  these  surfaces  do  not  have  perfect  form,  we  numerically 
fit  a perfect-form  feature. 

ISO  TC  213  defines  several  types  of  features.  For  example  a truncated  portion  of  a 
cylindrical  surface  on  a nominal  geometrical  model  of  a part  is  called  nominal  integral  feature 
(see  Figure  3)  . Its  axis  is  named  the  nominal  derived  feature  to  highlight  the  fact  that  the 
axis  is  derived  from  a cylindrical  surface  that  is  integral  to  the  part  boundary.  The  nominal 
geometrical  model  is  only  a concept  in  the  mind  of  a designer.  An  actual  realization  of  the 
part  will  have  a surface  that  only  approximately  corresponds  to  the  cylindrical  feature,  this  is 
called  the  real  integral  feature.  The  real  integral  feature  has  an  infinite  number  of  points  and 
for  the  purpose  of  measurements  only  a finite  subset  of  these  are  used.  These  finite  number  of 
points  sampled  on  the  real  integral  surface  form  what  is  called  the  extracted  integral  feature. 
Then  a perfect  form  surface,  such  as  a full  cylinder,  is  fitted  to  the  sampled  points  to  form 
what  is  termed  the  associated  integral  feature.  In  ISO/TC213,  the  term  association  is  used 
synonymously  with  fitting.  Finally,  associated  derived  feature  is  the  axis  or  other  appropriate 
entity  that  is  derived  from  the  fitted  surface. 

Deriving  a Distribution 

As  mentioned  already  we  fit  a perfect-form  feature  to  an  actual  feature.  As  shown  in  Figure  4, 
two  types  of  fits  are  popular:  Gaussian  and  Chebychev  fits. 

A Gaussian  fit  minimizes  the  squares  of  deviations  of  the  sensed  points  from  the  substitute 
element  and  the  solution  is  unique.  The  Chebychev  method,  on  the  other  hand,  uses  a min- 
max  criteria.  It  minimizes  the  maximum  deviation  of  the  sensed  points  from  the  substitute 
element.  The  choice  of  fit  may  be  dictated  by  the  function  of  the  part.  For  example,  if 
assemblability  is  the  major  functional  criteria,  then  Chebychev  fit  is  preferred.  The  planned 
ISO  10360  [27]  uses  Gaussian  fit  because  it  gives  a unique  solution. 
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Figure  3:  Various  types  of  features  defined  by  ISO  TC  213  ([71]) 


Acceptable  Distributions 

In  the  previous  section  we  discussed  the  process  of  mapping  a population  of  actual  features 
to  the  distribution  of  one  or  more  random  variables  of  interest.  In  this  section  we  discuss  the 
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Figure  4:  Extracting  parameters  or  actual  values  by  fitting  ([72,  73]) 


acceptable  class  of  population  of  actual  features  (parts),  by  defining  which  statistical  distri- 
butions are  acceptable.  In  [72,73],  three  ways  of  specifying  acceptable  class  of  distributions 
are  proposed  by  defining  statistical  tolerance  zones.  We  briefly  discuss  them  here.  For  more 
details,  we  refer  the  interested  reader  to  [72,  73]. 

2.5.2  Statistical  Tolerance  Zones 

The  statistical  tolerance  zones  are  of  two  types:  (1)  parametric  statistics,  and  (2)  nonpara- 
metric  statistics.  In  parametric  statistics,  the  population  parameter  zones  are  defined  in  the 
space  of  parameters  of  the  population.  Currently,  only  the  first  two  (central)  moments  of  the 
distribution  are  considered  in  the  population  parameter  space.  In  the  nonparametric  statis- 
tics, the  distribution  function  zones  are  defined  using  the  whole  distribution.  Conceptually, 
nonparametric  methods  are  more  general  than  parametric  methods. 

Parametric  Statistics 

Using  the  process  capability  indices  explained  in  Appendix  A (page:  67),  a process  is  sta- 
tistically controlled  by  specifying  lower  bounds  for  Cp  and  CPk,  and  an  upper  bound  for  the 
centering  index  (k  = Cc,  as  defined  in  [72]).  Essentially  we  need  to  specify  three  number 
P,  K and  F , so  that 


Cp  > P (lowerbound) 


(2) 
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Cpk  > K 
Cc  < F 


(lowerbound) 

(upperbound) 


(3) 

(4) 


The  condition  4 implies  Cpk  > (1  — F)Cp.  The  following  zones  fall  under  this  type: 

1.  Cp  Cpk  zone, 

2.  fj,  — a zone,  and, 

3.  Cpm  zone. 

Zone  in  Cpk  — Cp  Plane 

Each  value  of  the  triplet  (P,  K,  F)  corresponds  to  a tolerance  zone  in  the  Cpk  — Cp  plane. 
We  are  interested  in  the  first  quadrant  of  Cpk  — Cp  zone  where  both  Cpk  and  Cp  are  positive 
and  0 < F < 1.  The  tolerance  zone  is  the  intersection  of  the  four  half-planes 


cp 

> 

p 

Cpk 

> 

K 

cpk  - cp 

< 

0 

cpk 

> 

(1  - F)CP 

(5) 

Figure  5 shows  the  tolerance  zone  for  the  general  case  of  (P.  K , P)  and  the  Motorola-six 
sigma  case  (P  = 2,  K — 1.5,  P = 0.25).  In  case  we  do  not  have  lower  bound  for  Cpk , we  can 
use  the  condition  given  by  the  equation  (5)  above.  In  such  cases,  the  shape  of  the  zone  is 
similar  to  the  six-sigma  case. 


Cp— Cpk  Zone:  General  case  Cp— Cpk  Zone:  Motorola  Six— sigma  case 


Figure  5:  Tolerance  zone  in  Cpk  — Cp  plane 
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Zone  in  /i  — cr  Plane  Using  Cpk,  Cp 

The  triplet  (P,  K , F ) gives  rise  to  a different  shape  to  the  above  tolerance  zone  in  the  Qu  - a) 
plane.  Since  a > 0,  this  tolerance  zone  is  the  intersection  in  the  upper  half-plane  of  the 
following  half-planes. 

ri  USL-LSL  \ d ^ USL-LSL 

- 6o-  - r ^ ° - 6 P 

r — min{USL-y.,ti-LSL}  ^ jy-  . r-LSL  j - n+USL 

Upk  — 3a  — ^ a — 3/^  dna  a — 3K 

Cc  = < p =>  MPP  — Fd  < fj,  < MSL  + Fd 

where  M5  = \{LSL  + USL)  and  d — \{JJSL  — LSL).  Figure  6 illustrates  both  the  general 
and  the  six-sigma  case.  For  the  same  triplet  P,  iF,  P,  each  point  in  the  \i  — a tolerance  zone 
maps  to  a point  in  the  tolerance  zone  in  the  Cpk  — Cp  plane;  and  each  point  in  the  CPk  — Cp 
tolerance  zone  is  the  image  of  two  points  in  the  corresponding  /r  — a tolerance  zone  (from 
the  definition  of  CPk). 


Zone  in  n — a Plane  Using  Cpm 

We  use  the  index  Cpm  to  define  the  tolerance  zone  in  ^ — cr  plane  as  follows: 

Cpm  = L’SL  ~ LSL  = >M^{a2  +.{fi  - Tf)  < R 2 
6 ^ + (n  - Tf 

where  R2  — iUS3^2L'1  • We  get  the  tolerance  zone  as  a semi-circle  as  shown  in  the  Figure  7. 
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Figure  7:  Tolerance  zone  in  /r  — a plane  using  Cpm 

Nonparametric  Statistics 

In  the  previous  section,  we  discussed  tolerance  zones  that  are  defined  in  the  space  of  pa- 
rameters of  the  population,  namely  the  first  two  (central)  moments  of  the  distribution.  In 
the  nonparametric  statistics,  the  distribution  function  zones  are  defined  using  the  whole 
distribution.  If  the  measured  characteristic  X is  normally  distributed  and  if  we  plot  the  cu- 
mulative distribution  function  (CDF)  of  X on  normal  probability  paper,  it  will  be  a straight 
line.  The  slope  of  the  line  is  l/a  and  the  ^-coordinate  of  the  intersection  of  this  line  with 
y = 0.5  line  is  the  mean  (for  normal  distribution  mean  is  equal  to  median).  Each  line  with 
positive  slope  corresponds  to  a unique  normal  distribution. 

Each  point  in  the  /r  — a plane  corresponds  to  a line  in  the  normal  probability  plot  and  a 
line  segment  in  the  /i  — a plane  corresponds  to  a zone  in  the  normal  probability  plot  obtained 
by  sweeping  a line.  If  the  line  segment  in  /r  — a has  a non-zero  slope,  then  this  sweep  is  a 
pure  rotational  sweep  of  a line  in  the  normal  probability  paper.  The  horizontal  line  segment 
in  — a translates  to  a line  in  normal  probability  plot.  For  complete  details  refer  to  [72]. 
The  collection  of  normal  CDFs  which  belong  to  such  a tolerance  zone  is  independent  of  the 
choice  of  plotting  in  a normal  probability  paper.  This  collection  is  called  the  CDF  tolerance 
zone  [72].  Figure  8 shows  a typical  CDF  tolerance  zone. 

Composition  of  Statistical  Tolerance  Zones 

Until  now  we  have  discussed  statistical  tolerance  zones  for  piece  parts.  Now,  we  extend  this 
to  an  assembly  level  tolerancing.  In  order  to  compute  the  tolerance  zones  for  Y we  need  to 
formulate  the  rules  of  composition  of  statistical  tolerance  zones,  i.e.,  we  need  mathematical 
rules  for  combining  part  level  variations  to  produce  assembly  level  variations.  Specifying  in- 
dividual part  tolerances  statistically  as  axis-parallel  rectangular  zones  in  the  (/r,  a)  plane  will 
result  in  simple  evaluation  of  the  mathematical  rule.  In  the  general  case,  such  as  triangular 
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CDF  Tolerance  Zone 


17.55  17.56  17.57  17.58  17.59  17.6  17.61  17.62  17.63  17.64  17.65 


Figure  8:  CDF  Tolerance  Zone 


statistical  tolerance  zones,  the  rule  of  composition  involves  computing  Minkowski  sum  [74]  of 
objects  bounded  by  curve  segments  even  under  the  assumption  of  statistical  independence. 
The  rule  of  composition  for  CDF  tolerance  zones  is  totally  different  and  involves  advanced 
mathematical  computations.  The  working  group  WG  13  of  ISO  TC213  is  evaluating  this 
approach. 

Let  us  discuss  the  composition  of  statistical  tolerance  zones  for  a linear  case  of  ARF 
Y — a0  + a iX\  + a2X 2 + • • ■ + anXn,  where  a\s  are  some  constants.  We  also  assume  that 


X[s  are  independent  so  that 


= E a2ia 


(6) 


where  (/Ti,cr;)  are  the  mean  and  standard  deviation  of  Xx  respectively. 

Our  task  is  to  compute  the  statistical  tolerance  zone  for  Y.  To  do  this,  in  [74],  they  have 
defined  statistical  tolerance  zone  in  (/i  - a2)  plane  and  use  (6)  to  reduce  the  composition  of 
statistical  tolerance  zone  s to  simple  Minkowski  sum.  Refer  [74]  for  more  details. 


Statistical  Tolerance  Zones:  Summary 

Given  an  arbitrary  normal  distribution  for  the  variable  X,  we  can  check  if  the  population 
is  acceptable  by  checking  any  of  the  following:  (i)  if  the  tuple  (Cp Cp)  is  within  CPk  — Cp 
zone,  or  (ii)  if  an  ordered  pair  (/i.  — a)  is  within  the  (i  — a zone,  o (iii)  if  the  CDF  of  X lies 
within  the  CDF  zone.  If  X ~ a)  then  all  the  above  discussed  zones  are  equivalent;  that 
is,  the  normal  population  is  acceptable  according  to  one  interpretation,  if  and  only  if  it  is 
acceptable  in  the  other  two  interpretations  also  [72].  In  order  to  apply  this  methodology  for 
general  cases,  like  non-normal  and  unsymmetric,  we  need  to  generalize  the  above  tolerance 
zone  concepts.  This  may  be  necessary  for  the  interpretation  of  statistical  tolerancing  when 
applied  to  geometric  tolerances,  such  as  perpendicularity.  In  this  case  we  would  ideally  like 


32 


to  produce  cylindrical  features  with  zero  perpendicularity  tolerance.  This  implies  that  the 
target  should  be  LSL.  We  can  enrich  the  class  of  distributions  by  dropping  the  normality 
assumption.  This  may  be  necessary,  because  perfect  normality  is  seldom  realized  in  practice. 
Even  in  theory  the  normal  distribution  might  not  be  expected  for  certain  tolerances  (for 
example,  perpendicularity). 

If  the  process  is  centered  and  the  target  is  the  mid-point  of  the  specification  limits,  then 
the  indices  are  equivalent.  If  the  last  two  assumptions  are  relaxed  (that  is  unsymmetric  and 
non-normal)  care  should  be  taken  in  interpreting  the  values  of  the  indices.  But  as  long  as 
the  distributions  are  determined  by  their  first  two  moments  these  indices  will  serve  as  useful 
tool  to  assess  the  process  performance.  Figure  9 summarizes  the  various  zones. 


Parameter(s) 


Tolerance  Zone 


Parametric  Statistics 


Cpl 


Remarks 


Cpk-Cp 

zone 


More  appealing  to  process  and  quality  engineering. 
We  need  to  specify  (P,K,F) 


Equivalent  to  the  above.  Easier  to  perform 
statistical  tolerance  analysis  and  synthesis. 

In  case  1,  we  need  to  specify  (P.K.F)  and  in 
case  2 we  need  to  specify  lower  bound  for  Cpm 


2 

ixa 


zone 


Under  the  assumption  that  X’s  are  independent 
the  composition  of  ST  Zones  is  given  by  simple 
Minkowski  sum. 


Nonparametric  Statistics 


Pr(x  <=  X) 


CDF  zone 


X 


Greater  flexibility  and  generality  and  can 
accomodate  non-normal  and  unsymmetric 
situations.  If  the  X - N(p.,cy  ) and  process  is 
statistical  control,  then  all  the  zones  are 
equivalent. 

For  non-normal  cases,  as  long  as  the  distributions 
are  determined  by  the  first  two  (cenral)  moments, 
first  two  zones  agree  exept  for  minor  exeptions. 


Figure  9:  Various  Tolerance  Zones:  Summary 
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3 Function- Assembly-Behavior  Model 

In  this  report  our  main  aim  is  to  develop  an  integrated  comprehensive  function,  assembly 
(artifact)  and  behavior  model.  To  do  this,  we  need  to  develop  appropriate  assembly  models 
and  tolerance  models  that  would  enable  the  designer  to  incrementally  understand  the  build- 
up or  propagation  of  tolerances  and  optimize  the  layout,  features,  or  assembly  realizations. 
This  will  ensure  ease  of  tolerance  delivery.  In  [53],  a multi-level  approach  called  Design  for 
Tolerance  [DFT]  process  was  proposed.  The  DFT  enables  tolerancing  to  be  addressed  at 
successive  stages  of  design  in  an  incremental  fashion.  In  Section  4,  we  address  the  effective 
use  of  the  FAB  and  DFT  model  for  design  tolerancing  starting  from  the  conceptual  stage 
of  the  design  and  evolving  continuously  throughout  the  entire  design  process  to  the  final 
detailed  design. 

3.1  Literature  Review 

Attempts  to  describe  an  assembly  model  in  a CAD  environment  have  been  made  by  a number 
of  researchers  [21, 41, 62,  67,  70]  in  order  to  carry  out  assembly  analysis,  kinematics  analysis, 
and  tolerance  analysis.  Five  major  issues  related  to  automated  assembly  have  been  studied 
elaborately.  These  are:  1)  assembly,  component  and  feature  geometry,  and  topology  repre- 
sentation; 2)  identification  and  establishment  of  functional  relationships  between  components 
of  an  assembly;  3)  identification  of  precedence  relationships;  4)  assembly  sequence/plan  gen- 
eration; and  5)  assembly  analysis.  Since  mating  relations  establish  geometric  constraints 
between  components,  several  researchers  have  focussed  on  characterizing  mating  conditions 
and  the  impact  on  tolerance  synthesis  and  analysis.  Lee  and  Gossard  [47]  suggested  the 
inclusion  of  mating  information  in  a data  structure  representing  an  assembly  for  CAD/CAM 
systems.  Work  by  Roy  et  al.  [62, 67]  suggests  that  an  assembly  can  be  decomposed  into  com- 
ponent features  eventually;  thus  mating  conditions  among  components  can  be  represented  as 
a relation  among  features  by  a ‘modified  functional  relationship  graph’.  Srikanth  et  al.  [70] 
investigated  the  mating  relations  among  parts  and  developed  a ‘relation’  operator  to  repre- 
sent mating  conditions  in  Relational  Graph,  which  is  a graph  for  representing  an  assembly. 
This  representation  can  be  applied  for  dimensional  tolerance  analysis  of  an  assembly,  but  not 
for  geometrical  tolerances  among  features  of  the  component  parts.  Raj  an  et  al.  [61]  report 
the  development  of  assembly  representations  that  capture  mating  constraint  information 
generated  during  the  detailed  design  stage.  Important  mating  constraints,  joint  attributes, 
and  a relational  structure  are  explicitly  embedded  within  the  hierarchical  assembly  structure 
to  facilitate  assembly  analysis.  The  joint  kinematics  information  and  component  degrees  of 
freedom  available  from  the  assembly  model  can  be  used  in  formulating  a mathematical  model 
for  tolerance  analysis  and  synthesis. 

Wang  and  Ozsoy  [81]  have  described  representation  of  assemblies  for  automatic  gener- 
ation of  tolerance  chains.  The  assembly  is  represented  in  an  elaborate  data  structure  with 
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information  on  assembly  decomposition;  (4x4)  transformation  matrix  for  each  instance  of  a 
component/sub-assembly  mating  feature;  mating  conditions  (against,  parallel,  fit),  dimen- 
sions and  tolerances  of  mating  features,  etc.  A tolerance  chain  for  any  given  assembly  can  be 
generated  algorithmically  using  the  above  information.  Tolerance  analysis  uses  these  chains. 
The  representation  does  not  need  geometric  data  but  cannot  be  used  in  early  stages  of  design 
due  to  the  nature  of  information  required  to  complete  the  data  structure. 

In  [80],  a graph-based  representation  scheme-called  the  Tolerance  Network  (TN)-  was 
proposed  to  represent  the  geometric  dimensioning  and  tolerancing  (GD&T)  specifications  in 
a feature-based  product  model.  Tolerances  are  treated  as  variational  geometric  constraints 
and  the  authors  explain  how  the  tolerance  network  (symbolic  constraint  graph)  represents 
the  designer's  intentions  regarding  functional,  manufacturing,  and  quality  control  of  the 
product.  In  their  paper  the  authors  have  considered  ‘assembly  graph’  (an  undirected  graph 
with  parts  as  nodes  and  connection  constraints  as  arcs)  as  an  assembly  model  for  tolerancing. 
The  TN  can  be  a good  candidate  for  tolerance  representation. 

The  need  for  assembly  model  standards  is  one  of  the  major  suggestions  from  the  partici- 
pants of  the  national  workshop  at  NIST  [76].  The  Japanese  National  Committee  (JNC)  has 
proposed  a neutral  assembly  model  for  inclusion  in  STEP  [77].  An  extensive  review  of  as- 
sembly requirements,  representations,  and  methodologies  that  support  assembly  design  and 
planning  can  be  found  in[61] . Pratt  [51]  has  developed  a more  generic  geometric  constraint 
representation  for  use  in  describing  component  as  well  as  assembly  constraints. 

With  tolerance  analysis  as  the  main  objective,  Whitney,  Gilbert,  and  Jastrzebski  [82;  pro- 
posed a model  of  assembly.  This  model  contains  the  following  information:  mating  features 
that  build  up  the  assembly;  a graph  representing  the  mating  part;  an  underlying  coordi- 
nate structure  of  the  assembly;  and  homogeneous  (4x4)  matrix  transformations  to  represent 
dimensions  and  tolerances  of  each  part  (in  accordance  with  Y14.5M-1982  geometrical  toler- 
ancing standard).  The  transformations  represent  both  nominal  relations  between  parts  and 
variations  caused  by  geometric  deviations  allowed  by  tolerances.  These  transformations  can 
be  used  to  propagate  tolerances  through  an  assembly  so  as  to  compute  the  location  of  any 
designated  part  starting  from  a reference  part,  taking  into  account  variations  in  the  loca- 
tions, sizes,  and  shapes.  The  above  representation  can  potentially  be  used  in  the  early  design 
stage.  The  concept  of  Datum  Flow  Chain  (DFC)  advocated  by  Mantripragada  and  Whitney 
[49]  can  be  used  to  provide  information  necessary  for  locating  parts  accurately  with  respect 
to  each  other,  to  relate  the  datum  logic  explicitly  to  product  key  characteristics,  assembly 
sequences,  and  choice  of  mating  features  and  to  provide  information  needed  for  tolerance 
analysis. 

Ge  Qu  [60]  reported  a method  on  automatic  tolerance  allocation  (in  a linear  tolerance 
stack-up)  with  process  selection  based  on  hierarchical  assembly  model.  In  his  work,  an 
assembly  is  modeled  by  hierarchical  data  structure  and  assembly  relationships  between  com- 
ponents are  specified  by  mating  relations  through  the  mating  surfaces  of  each  component. 
Based  on  these  mating  relationships,  the  method  can  generate  dimension  graphs  (to  rep- 
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resent  dimension  vector  loops  in  the  assembly)  where  nodes  correspond  to  mating  surfaces 
and  edges  correspond  to  dimensions  (and  tolerances)  between  mating  surfaces.  The  dimen- 
sion graphs  provide  the  basic  mathematical  models  for  solving  tolerance  allocation  problems. 
This  research  work  addresses  linear  tolerance  stack-up  issues  for  only  three  types  of  fit  rela- 
tions: length  fits,  cylindrical  fits  and  housing-shafts.  It  uses  a tree  search-based  technique 
for  solving  discrete  tolerance  and  process  selection  problems  in  length  fits  and  a rule-based 
method  to  specify  tolerances  for  cylindrical  fits. 

What  is  missing  from  previous  works  is  a general  and  unified  representation  scheme  for 
geometric  tolerances  and  mating  conditions  across  parts  and  assemblies.  In  addition,  very 
little  consideration  has  been  given  to  the  significance  of  an  assembly  data  model  for  pro-active 
analysis  of  the  product  (especially  for  tolerance  synthesis  and  analysis)  as  design  progresses. 
In  this  paper,  we  discuss  the  integration  of  function,  assembly,  and  behavior  representation 
into  a comprehensive  information  model  (FAB  models).  To  do  this,  we  need  to  develop 
appropriate  assembly  models  and  tolerance  models  that  enable  the  designer  to  incrementally 
understand  the  build-up  or  propagation  of  tolerances  and  optimize  the  layout,  features,  or 
assembly  realizations.  This  will  ensure  ease  of  tolerance  delivery. 

3.2  Information  Requirements  for  Synthesis  and  Analysis 
of  Tolerances 

Significant  productivity  gains  can  be  achieved  if  we  consider  tolerance  information  through 
out  the  entire  product  design  cycle.  We  need  to  develop  appropriate  assembly  models  and 
tolerance  models  that  would  enable  the  designer  to  incrementally  understand  the  build-up 
or  propagation  of  tolerances.  We  also  need  to  optimize  the  layout,  features,  or  assembly 
realizations  so  as  to  ensure  ease  of  tolerance  delivery. 

Any  valid  assembly  sequence  has  a direct  influence  on  the  synthesis  of  tolerances.  At 
the  same  time,  it  should  be  kept  in  mind  that  any  tolerances  specified  also  directly  influence 
the  manufacturing  and  assembly  tasks.  The  tolerance  synthesis,  analysis,  and  assembly 
planning  tasks  cannot  be  treated  separately;  they  are  dependent  on  each  other.  Therefore, 
we  need  to  derive  some  kind  of  compromise  between  the  tolerance  design  and  the  various 
planning  tasks.  This  compromise  should  not  violate  the  established  rules  and  procedures  in 
assembly  planning.  It  should  be  able  to  take  account  of  the  lack  of  tolerance  information  in 
the  conceptual  stage  of  assembly  design  and  in  the  formulation  of  various  plans.  Different 
types  of  information  are  needed  at  different  abstract  levels  of  assembly  and  tolerancing.  At 
the  conceptual  level,  with  a view  towards  computerized  tolerance  synthesis  and  analysis,  we 
need  the  following  types  of  information:  a description  and  quantification  of  the  functional 
requirements  of  the  product-the  design  intent;  information  that  leads  to  the  selection  of 
assembly  configuration  (layout);  identification  of  datums,  features  and  geometric  relations 
between  features  that  will  influence  functional  and  assembly  requirements;  and  an  assembly 
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graph  that  provides  information  for  the  selection  of  location  logic  (in  deciding  the  manner 
in  which  parts  are  located  with  respect  to  one  another). 

During  tolerance  synthesis  and  analysis,  an  assembly  model  should  provide  information 
regarding  representation  of  assembly  design  in  terms  of  attributes  such  as  design  function, 
design  parameters,  geometric  characteristics,  manufacturing  features  and  processes,  and  as- 
sembly conditions  [63].  The  decomposition  of  overall  part  functional  requirements  (PFRs) 
and  their  mapping  into  different  product  specifications  are  important  in  the  context  of  prod- 
uct/assembly design  particularly  during  synthesis  of  geometry  and  tolerances.  The  most 
important  information  that  is  essential  for  tolerance  synthesis  is  the  geometric  characteris- 
tics of  each  component  and  inter-relationships  between  different  components  in  an  assembly. 
The  geometric  characteristics  of  a component  consist  of  size  (geometric  dimensions),  form 
(type  of  geometry,  i.e.  planar  or  cylindrical  surfaces)  and  profile  (nature  of  surface)  while 
inter-relationships  consist  of  relative  positioning  of  components. 

For  tolerance  analysis  in  the  detailed  engineering  stage,  an  assembly  model  should  provide 
kinematics  inputs  for  assembled  components  that  are  small  variations  of  component  dimen- 
sions around  their  nominal  values,  and  detailed  constraint  information  on  the  interaction 
between  mating  parts.  These  constraints  also  serve  as  functions- assembly  response  function 
(ARF)-by  which  assembly  variation  may  be  studied.  At  this  stage,  the  assembly  model  pro- 
vides information  on  dimension  vectors  that  is  used  to  represent  various  dimensions.  Each 
vector  represents  either  a component  dimension  or  kinematically  variable  dimension  that 
can  be  arranged  in  chains  or  loops  and  represents  those  dimensions  that  stack  together  to 
determine  resulting  assembly  dimensions.  The  other  model  elements  include  datum  refer- 
ence frames,  assembly  and  component  tolerance  specifications,  feature  datums  and  geometric 
feature  tolerances,  and  kinematics  joints  that  describe  motion  constraints  at  the  points  of 
contact  between  mating  parts.  It  should  also  provide  models  of  assembly  variations  that 
include  1)  dimensional  variations,  2)  geometric  variations,  and  3)  kinematics  variations. 

The  assembly  model  introduces  more  than  just  assembly  (and  manufacturing)  informa- 
tion for  tolerance  analysis  and  synthesis.  It  serves  other  functions,  such  as  the  functional 
analysis  of  the  whole  assembly,  and  assembly  sequence  issues.  The  assembly  model  provides 
detailed  knowledge/information  of  part  geometry  and  function. 

3.3  Information  Models  for  the  Design  for  Tolerance 

In  this  section,  we  will  first  describe  the  representational  issues  of  functional  models  and  the 
representation  of  the  FAB  model.  Then,  we  will  discuss  part  models  in  the  assembly  and 
applicable  tolerance  models  for  both  the  part  and  assembly. 

In  the  early  design  phase,  most  of  the  design  decisions  taken  are  concerned  with  the 
desired  characteristics  and  the  overall  functions  of  the  assembly.  By  function  we  mean  an 
abstract  formulation  (or  definition)  of  a task  that  is  independent  of  any  particular  solution. 
In  this  phase,  the  abstract  functional  specification  of  an  assembly  is  transformed  into  a 
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physical  description.  In  the  later  design  phases,  the  physical  decisions  that  are  made  in  the 
earlier  phases  are  elaborated  to  ensure  that  they  satisfy  the  specified  functional  requirements 
and  life  cycle  evaluation  criteria.  In  order  to  manipulate  the  function  information,  a func- 
tional data  model , that  describes  the  functional  information  throughout  the  design  cycle,  is 
needed.  Appropriate  reasoning  modules  can  then  interoperate  and  extract  functional  infor- 
mation during  the.  decision-making  processes  from  the  functional  data  model  (as  geometric 
reasoning  modules  query  data  from  the  product  data  model , that  is  the  CAD  model,  for  shape 
information). 

3.3.1  Function  Representation  and  Data  Model 

Depending  on  the  design  phase,  different  types  of  functions  exist  at  different  levels  of  ab- 
straction. In  the  conceptual  design  phase,  functions  are  usually  independent  of  the  product 
domain,  whereas  in  later  design  phases,  when  the  functions  are  detailed,  they  become  in- 
creasingly dependent  on  the  product  domain  that  has  been  selected.  Therefore,  we  need 
to  adopt  an  appropriate  formal  function  representation  scheme  that  will  be  helpful  when 
modeling  the  overall  function  of  the  assembly  in  the  conceptual  design  stage,  and  will  also 
be  useful  to  model  a smaller  sub-assembly,  particularly  as  component  and  feature  levels  are 
approached.  The  representation  should  also  have  unambiguous  semantics  in  order  to  perform 
analysis  on  (i.e.  to  manipulate)  the  functional  description.  Many  researchers  [5, 12,44]  have 
already  investigated  different  types  of  functional  representation  schemes.  Among  these,  two 
types  of  representation  are  particularly  relevant  in  the  domain  of  electro-mechanical  prod- 
ucts: 1)  model  of  the  flows  between  inputs  and  outputs  of  a product,  and  2)  syntactic 
languages,  which  use  a grammar  consisting  of  verbs  and  nouns  to  describe  a product.  Refer 
to  [5, 12,35,44,49]  for  a detailed  discussion  on  these  techniques.  Both  representations  have 
limited  use  in  the  synthesis  of  mechanical  devices.  They  are  not  able  to  support  applications 
over  the  complete  product  life  cycle.  Our  functional  data  model  can  represent  the  functional 
description  of  a product  through  out  the  entire  product  design  cycle. 

For  computer  representation  of  function  and  behavior,  we  adopt  the  object-oriented  de- 
sign modeling  language  that  has  been  developed  as  part  of  the  Design  Repository  Project 
at  the  National  Institute  of  Standards  and  Technology  [52].  A Function  has  a set  of  inputs 
and  outputs,  and  may  have  other  attributes  or  relationships  that  specify  sub-functions  and 
composite  functions.  Functions  are  treated  as  the  instantiations  of  classes  of  functions.  The 
Function  hierarchy  consists  of  four  subclasses  of  function:  Transform_function,  Con- 
vey .function,  Supply_function  and  ControlJuncion  (for  a detailed  description,  refer 
[52]).  Inputs  and  Outputs  of  functions  are  described  in  terms  of  sets  of  instances  of  a sub- 
class of  the  general  Fluent  class.  Fluents  include  both  physical  and  abstract  phenomena 
that  are  associated  with  inputs  and  outputs  to  functions,  such  as  motion,  force,  heat,  liquids, 
current,  and  so  on.  The  Fluent  class  is  divided  into  two  main  subclasses:  Abstract  .fluent 
and  PhysicalJluent.  The  Abstract  .fluent  class  specifies  things  like  motion  and  electric- 
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ity  while  Physical -fluent  class  specify  things  like  gases  and  liquids.  All  Fluent  classes 
have  slots  for  explicit  references,  to  which  artifact  the  fluent  is  flowing  from,  and  where  it 
is  flowing  to.  The  set  of  subclasses  for  the  general  Function  and  Fluent  classes  together 
describe  the  semantics  of  the  functional  and  behavioral  aspect  of  the  object-oriented  design 
modeling  language. 

Function- Assembly-Behavior  Model 

Figure  10  depicts  a typical  FAB  data  model  for  capturing  the  product  development-related 
issues  from  the  conceptual  design  stage  to  the  detailed  assembly  building  process.  The  main 
purpose  of  the  proposed  aggregate  structure  of  the  function,  behavior  and  assembly  together 
in  this  data  model  is  to  support  conceptual  design  as  well  as  design  for  manufacturing  and 
assembly,  starting  from  an  early  design  stage.  It  starts  with  defining  the  design  context 
that  represents  a particular  design  domain  and  consists  of  the  design  tasks  (i.e. , the  Goal) 
relevant  to  the  current  design  assignment,  the  decisions  that  have  been  made,  the  design  re- 
quirement description  and  the  artifact  that  has  been  developed  as  part  of  the  design  process. 
The  design  requirements  have  three  attributes:  description,  weight  and  dictatedJby. 
The  description  records  the  requirement  description  that  covers  input  specifications  by 
the  customer,  objectives,  and  specified  constraints.  Following  a formal  method  like  Quality 
Function  Deployment  (QFD),  customer  requirements,  stake-holder  requirements,  operational 
requirements,  and  other  statutory  requirements  are  developed.  These  requirements  can  be 
further  categorized  into:  1)  spoken  expected  requirements  (i.e.  directly  specified  by  the  cus- 
tomer); 2)  unspoken  expected  requirements  (i.e.  automatically  expected  by  the  customer); 
and  3)  unspoken  unexpected  requirements  (i.e.  no  direct  need  but  it  will  be  nice  to  have 
them).  The  weight  attribute  describes  whether  the  requirement  is  a demand  or  wish.  The 
product  requirements  are  finally  fulf  illed_by  a set  of  functions. 

The  Function  (a  function  means  what  the  system/product  is  for)  is  associated  with 
the  transformation  of  an  input  physical  entity  into  an  output  physical  entity  by  a part. 
A function  has  the  attributes:  name,  inputs,  outputs,  dictatedJby,  and  achievedJby. 
The  transformation  between  the  physical  entities  (or  fluents  in  inputs  and  outputs)  is  con- 
trolled by  the  behavior  of  the  Artifact  (i.e.,  assembly).  This,  in  turn  is  governed  by 
part /sub-assembly /assembly  geometry,  either  by  itself  or  in  association  with  other  parts/sub- 
assemblies in  an  assembly,  and  the  physical  laws  which  are  associated  with  those  physical 
entities.  The  fulf  illedJby  attribute  of  the  function  relates  it  with  a (or  a set  of)  product 
requirements.  The  fulf  illedJby  and  dictated_by  relationships  between  product  require- 
ments and  functions  thus  facilitate  both  top-down  and  bottom-up  of  design  synthesis  pro- 
cesses. An  Artifact  which  can  either  be  an  assembly,  sub-assembly,  individual  part,  feature, 
or  even  any  low-level  geometric  entity  (i.e.,  face,  edge,  or  vertex)  achieves  the  Function. 

The  assembly /part  behavior  is  achieved  through  certain  functional  relationships  associ- 
ated with  its  sub- assembly /features  [20].  In  the  conceptual  design  phase,  each  functional 
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Figure  10:  Function  Assembly  Behavior  Model 


relationship  is  translated  to  its  spatial  equivalents  using  appropriate  physical  laws.  A part 
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achieves  its  behavior  through  the  interactions  of  its  geometry  or  environment.  Thus,  the 
utility  of  part  behavior  can  only  be  realized  when  it  is  decomposed  into  an  equivalent  part- 
function  model.  The  design  solutions  are  obtained  as  a composition  of  physical  solutions 
and  technical  solutions.  A physical  solution  is  the  first  step  to  realize  behavioral  require- 
ments and  is  governed  by  physical  laws.  The  technical  solutions  are  technical  manifestations 
of  the  physical  solution.  A technical  equivalence  (i.e.  a physical  product)  determines  the 
final  product/artifact  and  defines  final  target  values  that  are  expected  to  achieve  the  desired 
function. 

The  assembly  data  model  contains  assembly-specific  information  of  a product  that  not 
only  helps  in  detailed  assembly  analysis  (i.e.,  determining  assembly  sequence,  etc.)  but  also  in 
synthesizing  the  system’s  design  and  studying  its  behavior.  Its  attributes  are:  key  assembly 
characteristics,  required  assembly  equipment  and  accessories,  sub-assembly/part, 
assembly  tolerances  and  assembly  kinematics  variations.  The  product’s  important 
geometric  features  and  material  property  information  (directly  derived  from  the  customer  re- 
quirements) that  are  highly  constrained  are  given  by  key  assembly  characteristics  [49]. 
These  assembly  key  characteristics  cannot  be  violated  in  any  situation  to  achieve  a desired 
level  of  the  product’s  performance,  function,  and  form.  The  description  and  specification 
of  required  assembly  equipment  and  accessories  are  necessary  for  the  assembly  analysis  at 
the  detailed  design  level.  The  sub- assembly /part  defines  the  composition  of  an  assembly. 
Assembly  tolerances  are  applied  to  the  assembly  dimensions  which  are  kinematics  assembly 
variables,  that  arise  during  the  assembly  of  components.  It  should  be  noted  that  while  com- 
ponent dimensions  are  subjected  to  random  process  variations,  the  assembly  dimensions  are 
not  manufacturing  process  variables;  these  are  kinematics  assembly  variables.  Kinematics 
variations  occur  at  assembly  time,  whenever  small  adjustments  between  mating  parts  are 
required  to  accommodate  dimensional  or  form  variations. 

The  sub-assembly /part  has  five  major  attributes:  name,  identification  number,  a set 
of  features,  relationship  with  other  sub-assembly  or  part  and  a pointer  to  its 
parts’  CAD  database  (that  supplies  geometrical,  topological,  variational  (tolerance)  and 
other  technical  information).  The  parts’  CAD  database- contains  detailed  constructive  solid 
geometry/boundary  representation  (CSG/B-Rep)  information.  A feature  is  composed  of 
feature  identification  number,  feature  datum,  feature  locations  and  orientation,  geometric 
feature  tolerance,  process  characteristics  for  manufacturing  the  feature,  and  its  spatial  rela- 
tion with  other  feature(s). 

Since  dimensional  variations  (lengths,  angles,  etc.),  geometric  form,  and  feature  variations 
(flatness,  roundness,  angularity,  etc.)  are  the  result  of  variations  in  the  manufacturing  pro- 
cesses, the  process  characteristics  hold  key  attributes.  The  spatial  relation  of  a feature  with 
other  feature(s)  describes  1)  datum  reference  frames  (DRFs),  2)  type  of  relationships  among 
features  (i.e.,  feature  interaction  list),  and  3)  and  any  other  feature  characteristics  (e.g.,  face 
pointers  to  the  respective  B-Reps  of  the  CAD  database).  The  spatial  relations  (with  other 
feature(s))  may  include  constraints  such  as  articulation  representation  for  electro-mechanical 
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Figure  11:  Function  Assembly  Behavior  Model: Assembly 


part  at  the  feature  level  (e.g.,  position  coincidence,  rotation  axis  coincidence,  parallel  nor- 
mals, center  coincidence) 

In  the  sub- assembly /part  level  the  relationship  of  a sub-assembly /part  with  another  sub- 
assembly/part has  two  main  attributes:  relationship  type  and  mating  characteristics.  There 
are  various  relation  types:  contact,  detachable  type  rigid  attachment,  constraint  etc. 
Mating  characteristics  are  described  by:  1)  mating  sub-assembly /part  identification  number, 
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2)  identification  numbers  of  mating  features,  and  3)  mating  conditions.  Each  face  contains 
the  following  information  about  face-type  (planar,  cylindrical,  etc.): 

• dimension, 

• size_tolerance, 

• f orm_tolerance, 

• normal  vector, 

• process  model,  force_d.o.f, 

• kinematic_d . o . f , 

• surf ace.roughness,  and 

• f ace_interaction_list,  refer  [7]  for  details. 

Mating  conditions  are  described  by  mating  type,  joint  characteristics,  spatial  constraints, 
rules,  and  attributes.  Mating  type  can  be  any  of  the  following:  fit  (i.e.,  sliding,  clearance, 
transition,  and  interference), fasteners,  gear,  bearing,  etc.  These  can  be  further  described  in 
terms  of  some  basic  relationships:  orientation,  primary  mating  relations  (no-contact, 
touch-contact  and  overlap-contact),  size  relation,  intersects,  between,  contains,  and 
connect. 

In  this  representation,  assembly  joints  are  defined  by  two  joint  attributes:  rigid  and 
movable  [61] . Rigid  joints  require  only  the  position  (location  and  orientation)  information 
to  be  specified,  whereas  movable  joints  require  the  specification  of  position  and  kinematics 
information.  Kinematics  information  specifies  the  type  of  joint  (revolute,  cylindrical,  pris- 
matic, spherical,  etc.).  Note  that  this  information  is  explicitly  embedded  in  the  assembly 
data  model  for  the  purpose  of  information  completeness.  This  information  is  directly  usable 
in  the  detailed  assembly  analysis  phase.  Knowing  the  contact  relationships  between  various 
parts  would  facilitate  the  determination  of  joint  characteristics  in  an  automated  manner. 

Similarly,  spatial  constraints  may  also  be  explicitly  specified  or  may  be  implicitly  specified 
by  the  adjacency  of  component  surfaces.  Depending  on  the  nature  and  degrees-of-freedom, 
spatial  constraints  can  be  of  different  types:  completely  unconstrained,  single  planar,  multi- 
ple planar,  planar  motion,  constrained  planar  motion,  cylindrical,  uni-cylindrical,  prismatic, 
uni-prismatic,  revolute,  and  completely  constrained  (refer  to  [61]  for  a detailed  discussion). 

Rules  and  attributes  represent  other  necessary  technical  information  about  the  mating 
condition.  The  technical  information  describes  properties  of  the  joint  connection  and  re- 
strictions that  are  applied  in  the  connection  between  parts/sub- assemblies.  For  example,  in 
a rivet  joint,  the  strength  of  the  rivet  joint  is  described  by  attributes  and  rules  of  the  joint 
connection. 
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Figure  12:  Function  Assembly  Behavior  Model  - Mating  Characteristics 

The  Function- Assembly-Behavior  data  model  described  so  far  provides  information  for 
the  analysis  of  product  life-cycle  process.  Not  all  the  information  is  available  in  the  early 
conceptual  design.  The  information  content  of  the  data  model  evolves  with  the  product 
design  process. 

The  complete  object  oriented  class  definitions  for  Function,  Artifact  and  Behavior  are 
under  development  and  will  be  reported  in  a future  NIST  publication. 
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Figure  13:  Function  Assembly  Behavior  Model-  Mating  Types 

4 FAB-DFT  Integration 

In  the  previous  section  we  discussed  and  explained  the  object  oriented  architecture  of  the 
FAB  model.  In  this  section  we  discuss  how  a designer  can  use  this  FAB  model  for  design 
tolerancing.  In  [53]  a dynamic  architecture,  called  Design  For  Tolerancing  (DFT),  has  been 
proposed  to  help  the  designers  to  make  tolerance-related  decisions  from  the  conceptual  design 
to  detailed  design.  The  basic  motivation  of  the  DFT  architecture  is  the  definition  of  a multi- 
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level  approach  that  enables  tolerancing  to  be  addressed  at  successive  stages  of  design  in  an 
incremental  fashion.  The  resulting  design  process  integrates  three  important  domains:  (1) 
design  activities  at  successive  stages  of  design;  (2)  assembly  models  that  evolve  continuously 
through  the  design  process;  and  (3)  methods  and  best  practices  for  tolerance  analysis  and 
synthesis. 

The  basic  motivation  for  defining  the  FAB  model  is  to  change  the  traditional  thinking  of 
design  as  form-to-function  transformation.  Present  CAD  systems  have  a wealth  of  tools  for 
generating  geometric  forms  for  the  objects.  This  forces  the  designers  to  decide  on  the  form 
first  then  think  about  the  function  later.  This  way  of  designing  has  the  potential  problem  of 
designing  sub-optimal  designs  that  are  a mere  modifications  of  existing  ones.  Present  CAD 
systems  cannot  help  the  designers  to  be  innovative.  Hence,  tools  for  helping  a designer  to 
think  in  terms  of  function  should  be  developed  and  form  should  result  from  function,  i.e., 
function-to-form. 

Knowledge-based  design  systems  implement  this  paradigm  by  first  focusing  on  the  sym- 
bolic aspects  of  design  and  later  mapping  the  symbolic  structure  to  a geometric  model.  They 
can  also  capture  the  various  semantic  relationships  between  design  objects.  Essentially, 
knowledge-based  systems  use  techniques  developed  by  artificial  intelligence  researchers  to 
capture  the  knowledge  of  expert  designers  in  a computer.  These  systems  help  the  engineer 
in  the  design  generation  process. 

The  proposed  FAB  model  goes  one  step  further  to  include  behavior  modeling  aspects 
into  function-to-form  method  of  design.  The  various  classes  defined  in  the  object  oriented 
model  of  the  FAB  will  help  a designer  to  study  design  activity,  providing  various  viewpoints. 
For  example,  a designer  can  study  assemblability  analysis,  ease  of  tolerance  delivery,  quality 
and  supply  chain  etc.  Our  main  focus  is  on  tolerancing.  The  three  threads  of  the  DFT 
architecture  are  explained  in  Figure  14,  source  [53]. 

The  FAB  model  evolves  the  function  and  form  (assembly)  for  the  entire  design  life  cycle. 
The  FAB-DFT  integration  is  better  explained  through  an  example.  Before  we  discuss  the 
example,  we  discuss  briefly  about  various  research  efforts  on  assembly  process  models. 

In  the  literature,  several  researchers  have  presented  their  viewpoint  of  what  the  assembly 
design  process  should  be.  In  [53]  a brief  outline  of  some  viewpoints  that  emphasize  toler- 
ancing is  provided  using  the  SIMA  (Systems  Integration  for  Manufacturing  Applications) 
reference  architecture  formulated  at  the  National  Institute  of  Standards  and  Technology 
[3],  which  provides  a generic  specification  of  design-related  activities  for  electro-mechanical 
products.  Figure  4 shows  the  various  design  stages  and  activities  in  the  SIMA  reference  archi- 
tecture. Stage  All  ( Plan  Products ) involves  developing  the  idea  for  the  assembly  depending 
on  market  needs  and  customer  requirements  and  characterizing  it  in  terms  of  function,  tar- 
get price  range,  and  relationship  to  existing  product  lines.  In  Stage  A12  ( Generate  Product 
Specifications) , an  engineering  specification  for  the  assembly  is  formulated.  This  involves 
mapping  the  customer  requirements  into  engineering  requirements  and  refining  these  in  con- 
sideration of  the  relevant  laws,  regulations,  patents,  and  product  standards,  etc.  In  Stage 
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ASSEMBLY  MODELS  AND  REPRESENTATION 


DESIGN  FOR  TOLERANCE  PROCESS 


METHODS  AND  BEST  PRACTICES  FOR  TOLERANCING 


Figure  14:  Assembly  models,  design  process  stages,  and  tolerancing  tools 
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A1  : DESIGN  PRODUCT 


All  : Plan  Products 

A12  : Generate  Product  Specifications 

A13  : Perform  Preliminary  Design 

A131  : Develop  Functional  Decompositions 
A132  : Evaluate  and  Select  Decompositions 
A133  : Develop  Preliminary  Configurations 
A134  : Consolidate  Configurations 
A135  : Evaluate  Alternative  Designs 
A136  : Select  Design 
A14  : Produce  Detailed  Designs 
A141  : Design  System/Component 
A142  : Analyze  System/Component 
A143  : Evaluate  System/Component  Design 
A144  : Optimize  Designs 
A 145  : Produce  Assembly  Drawings 
A146  : Finalize  System/Component  Design 


Figure  15:  Design  stages  and  activities  in  the  SIMA  reference  architecture.  Source:  [3] 

A13  ( Perform  Preliminary  Design ),  the  assembly  design  problem  is  decomposed  into  a set  of 
component/subassembly  design  problems  and  specifications  are  developed  for  each  compo- 
nent/subassembly problem.  Interface  specifications  between  the  components/sub-assemblies 
are  developed  and  a preliminary  assembly  layout  is  created.  Finally,  in  Stage  A14  ( Produce 
Detailed  Designs ),  all  specifications  needed  to  completely  describe  each  subassembly  or  com- 
ponent are  produced.  This  includes  drawings  and  geometry,  materials,  finish  requirements, 
assembly  drawings,  and  fit  and  tolerance  requirements. 


4.1  Design  Tolerancing:  An  Incremental  Process 

Potentially,  tolerance  considerations  can  influence  the  decisions  taken  at  different  design 
stages,  in  increasing  level  of  detail.  Also,  the  decisions  taken  at  a particular  stage  influence 


48 


and  can  simplify  the  decisions  taken  in  the  downstream  stages.  Like  other  attributes  of 
a product  design,  tolerance  information  changes  over  time,  through  successive  stages  from 
product  planning  to  detailed  design  through  on-going  production.  Hence  a robust  tolerance 
representation  would  be  mutable  and  directly  related  to  the  evolving  product  design.  The 
incremental  refinement  of  design  tolerancing  processes  and  tolerance  representations  pro- 
ceeds in  symbiotic  fashion.  Consider,  for  example,  a tooling  design/build  process.  Both 
lead  time  and  cost  for  tooling  is  often  highly  dependent  on  the  tightness  of  the  tolerance 
requirements.  Scheduling  of  rough  cutting  for  a die  or  mold  can  typically  proceed  prior  to  a 
final  tolerance  specification,  but  the  finish  cut,  polishing,  etc.  must  proceed  afterward.  Con- 
versely, tolerance  specification  for  a snap-fit  in  a high-precision  injection-molded  part  must 
be  preceded  by  a decision  about  assembly  process  (e.g.,  manual  or  robotic).  For  complex  as- 
semblies with  many  parts,  the  timing  and  precedence  of  tolerancing  decisions  can  profoundly 
affect  scheduling  and  total  lead  time.  Analysis  and  synthesis  for  critical  tolerance  stack-ups 
is  clearly  related  to  process  plan  refinements.  There  are  opportunities  to  compress  cycle 
time  by  improved  modeling  prior  to  detailed  design,  but  compatible,  incrementally-refined 
representations  of  tolerances  and  manufacturing  processes  are  the  key. 

The  incremental  and  continuous,  ongoing  nature  of  tolerance  decision  making  enables  a 
natural  aggregation/decomposition  of  tolerancing  activities  as  the  design  matures.  Another 
way  of  viewing  this  is  in  terms  of  the  pruning  that  this  causes  at  successive  stages  in  the 
space  of  feasible  solutions  to  the  design  problem.  Early  on  in  the  design  process,  the  solution 
space  has  a staggering  cardinality  and  tolerancing  decisions,  if  taken  in  a continuous  ongoing 
fashion,  can  lead  to  substantial  early  reduction  in  the  space  of  possible  solutions  thus  making 
the  design  process  efficient.  Another  alternative  view  is  in  terms  of  marked  reduction  in 
design  iterations  or  design  rework.  In  this  sense,  design  for  tolerance  is  similar  in  spirit  to 
design  for  manufacturing/assembly  [9]  that  also  has  the  effect  of  dramatically  shrinking  the 
space  of  solutions  and  reducing  iterations.  Furthermore,  DFA,  DFM,  or  such  other  design 
related  strategies  may  have  close  coupling  with  tolerance-related  decisions  and  may  both 
influence  and  be  influenced  by  tolerancing  at  various  stages. 

4.2  Design  for  Tolerance:  A Multilevel  Approach 

The  first  two  stages  All  and  A12  of  the  SIMA  reference  architecture  and  also  the  early 
stages  of  other  assembly  design  processes  essentially  involve  mapping  customer  requirements 
into  product  ideas  and  product  specifications.  Tolerancing  is  not  directly  involved  in  these 
early  stages,  except  in  very  abstract  terms;  however,  these  stages  provide  critical  inputs  to 
the  tolerancing  decisions  in  the  rest  of  the  design  process. 

Thus  we  focus  on  Stage  A13  ( Perform  Preliminary  Design ) and  Stage  A14  ( Produce 
Detailed  Designs ) of  the  SIMA  reference  architecture.  We  divide  these  stages  into  the  fol- 
lowing four  tolerance-related  levels  (TR  Level)  and  develop  a four-level  approach  to  design 
tolerancing.  Note  the  difference  between  SIMA  stages  and  tolerance-related  levels  here. 
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• SIMA  Stage  A13:  Perform  Preliminary  Design 

- TR  Level  1:  Assembly  Layout  and  Configuration 

- TR  Level  2:  Location  Logic  and  Assembly  Features 

- TR  Level  3:  Assembly  Planning  and  Sequencing 

• SIMA  Stage  A14:  Produce  Detailed  Designs 

- TR  Level  4:  Detailed  Tolerance  Analysis  and  Synthesis 

These  levels  are  fairly  representative  and  generic  for  electro-mechanical  assemblies.  Neither 
the  number  of  levels  nor  the  description  of  the  individual  levels  is  to  be  viewed  as  being 
definitive.  Figure  16  captures  the  essence  of  this  architecture  for  DFT. 


4.2.1  TR  Level  1:  Assembly  Layout  and  Configuration 

Once  the  product  concept  is  known  and  engineering  specifications  are  generated  based  on 
the  key  characteristics,  TR  Level  1 of  the  proposed  process  can  commence.  TR  Level  1 in- 
volves decisions  regarding  the  preliminary  assembly  layout /configuration.  Such  decisions 
may  include:  rough  allocation  of  space,  number  of  sub- assemblies,  the  configuration  of  crit- 
ical sub-assemblies,  grouping  of  components  into  sub-assemblies,  and  rough  layout  of  the 
assembly. 

Using  the  FAB  model  we  need  to  generate  the  following  information  to  make  tolerance- 
related  decisions.  They  are: 

1.  rough  shapes/form  for  the  parts/sub-assemblies 

2.  parts  list 

3.  parts  location 

4.  layouts  and  configurations 

The  information  thus  generated  can  be  described  in  the  form  of  a liaison  diagram  (rela- 
tions between  parts  or  sub-assemblies),  a tree  (assembly  decomposition),  and  a partial  DFC 
(to  capture  whatever  location  logic  is  known  at  this  point).  Candidate  layouts  or  configura- 
tions can  be  identified  and  represented  using  these  models.  These  layouts  or  configurations 
and  related  manufacturing  process  selection  typically  might  differ  in  terms  of  ease  of  tol- 
erancing.  The  tolerancing  considerations  here  are  at  a coarse  level  and  may  be  directly 
influenced  by  customer  specifications.  To  effect  such  high  level  tolerancing  decisions,  aggre- 
gate level  manufacturing  process  capability  data  will  be  required  and  is  often  available  at 
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Assembly  Models  and  Tolerance  Information  DFT  Process 


Figure  16:  An  architecture  for  design  for  tolerance.  Legend:  LD  - Liaison  Diagram;  DFC 
- Datum  Flow  Chain:  HM  - Hierarchical  Model;  OOM  - Object  Oriented  Model;  AS  - Assembly 
Sequence;  PKCs  - Product  Key  Characteristics;  AKCs  - Assembly  Key  Characteristics;  ARF  - 
Assembly  Response  Function;  PCD  - Process  Capability  Data 
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this  point.  Simple  statistical  assumptions  and  probabilistic  calculations  can  be  used  at  this 
stage.  Also,  for  problems  such  as  manufacturing  process  selection,  the  formulations  such  as 
in  Section  2.2.3  can  be  used. 

4.2.2  TR  Level  2:  Location  Logic  and  Assembly  Features 

At  the  next  level  (TR  Level  2),  the  following  information  is  generated  using  the  FAB  model: 
assembly  response  functions  (approximate);  tolerance  requirements  at  interfaces  between 
major  sub-assemblies  and  parts.  The  relevant  process  capability  data  could  also  made  avail- 
able. The  decisions  here  are  concerned  with  the  location  logic  (how  to  locate  sub-assemblies 
and  components  with  respect  to  one  another)  and  with  choosing  the  appropriate  assembly 
features  to  go  with  the  location  logic.  The  choice  of  features  itself  might  depend  on  the 
assembly  sequence  (not  the  detailed  sequence  but  a precedence  specification  among  major 
assembly  steps).  The  DFC  model  is  suitable  to  capture  the  available/evolving  assembly 
information  here.  There  is  close  coupling  among  selection  of  features,  selection  of  assembly 
sequence,  and  creation  of  DFC.  Assembly  models  such  as  liaison  diagrams  are  also  relevant 
here.  If  the  assembly  is  of  Type  1,  then  the  assembly  features  are  predominantly  decided  by 
the  functional  requirements;  if  the  assembly  is  of  type  2,  then  the  choice  of  assembly  features 
is  an  important  problem  by  itself.  In  the  latter  case,  the  DFC  alone  will  not  be  adequate  to 
conduct  a tolerance  analysis.  A more  detailed  model  that  captures  the  tolerance  flow  at  this 
level  will  be  required.  Tolerance  analysis  here  can  tell  us  which  location  logic  is  better  from 
a tolerancing  viewpoint  and  which  set  of  assembly  features  would  best  accomplish  tolerance 
achievement.  This  stage  might  also  help  us  to  find  preliminary  target  values  and  tolerances 
for  individual  parts. 

Statistical  tolerance  analysis  methods  mentioned  in  the  preceding  sections  are  all  relevant 
here.  Determining  robust  nominal  values  and  preliminary  settings  of  tolerances  can  be 
accomplished  using  Taguchi  methods  or  Xerox  HPD  methodology. 

4.2.3  TR  Level  3:  Assembly  Planning  and  Sequencing 

We  proceed  next  to  TR  Level  3 where  the  detailed  assembly  response  function,  skeletal 
geometry  of  the  assembly,  assembly  features,  and  specification  of  parametric  or  geometric 
tolerances  of  individual  parts  and  features  are  again  generated  using  the  FAB  model.  The 
PCD  could  be  made  available  (if  applicable).  From  the  tolerance  specification,  one  may 
derive  (4  x 4)  matrix  transforms  for  the  nominals  and  variabilities  associated  with  the  parts 
[82].  The  decisions  here  could  be  with  respect  to  the  selection  of  the  detailed  assembly 
sequence  that  achieves  the  required  tolerance  specifications  in  the  best  possible  way.  The 
models  that  we  employed  in  the  previous  stage,  like  DFC  and  liaison  diagrams,  can  again 
be  used  here.  In  fact,  they  are  now  updated  with  richer  and  more  detailed  information. 
This  kind  of  representation  and  analysis  is  presented  in  [81],  where  several  data  structures 
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to  capture  tolerance-related  information  are  presented.  With  the  information  available  here, 
one  can  also  carry  out  tolerance  synthesis. 

4.2.4  TR  Level  4:  Detailed  Tolerance  Analysis  and  Synthesis 

TR  Level  4 corresponds  to  the  detailed  assembly  design  stage.  Here,  the  complete  assem- 
bly sequence  is  known;  geometric  data  about  the  parts  and  features  is  available;  detailed 
part  level  tolerance  requirements  are  known;  the  assembly  response  function  is  available  in 
complete  form;  and  low  level  process  capability  data  is  accessible.  Detailed  tolerance  anal- 
ysis and  synthesis  can  be  carried  out  here.  Most  tolerancing  studies  and  tolerancing  tools 
available  support  this  level  of  design. 

5 Current  and  Evolving  Standards:  Useful  Inputs  from 
FAB-DFT 

The  FAB  model  presented  in  this  report  can  eventually  lead  to  standards  for  tolerance 
and  assembly  representation  and  design  process  and  product  integration.  The  DFT  and 
FAB  integrated  model  will  help  the  designers  to  make  tolerance-related  decisions  from  the 
conceptual  design  to  detailed  design  and  throughout  the  entire  design  life  cycle. 

The  FAB  and  DFT  framework  uses  both  the  top-down  and  bottom-up  approach  to  the 
product  design.  The  hierarchical  assembly  model  defined  in  this  report  uses  features  as 
the  atomic  objects.  The  FAB  model  is  defined  in  such  a way  that  the  various  definitions 
and  concepts  that  are  already  in  use  (ISO  10303-47  [26],  ASME  Y 14.5  [1,2])  and  current 
evolving  standards  can  be  easily  incorporated  into  the  model.  The  FAB  model  can  help 
the  standards  community  in  specifying  the  representation  for  assembly  models,  tolerance 
representation,  and  assembly  level  tolerancing.  The  class  structure  (Function,  Behavior, 
Artifact)  defined  in  FAB  will  help  in  unifying  the  various  current  and  evolving  definitions, 
concepts,  and  standards.  In  the  following  sections  we  briefly  discuss  ISO  10303-47  and  the 
Geometric  Product  Specifications  (GPS)  standards  being  developed  by  ISO  TO  213. 

5.1  ISO  10303-47:  Shape  variation  tolerances 

This  part  of  ISO  10303,  [26],  supports  the  dimensioning  and  tolerancing  methods  defined  by 
the  ISO  1101  [31]  family  of  standards  for  the  dimensions  and  tolerances  on  the  engineering 
drawings.  These  methods  include  explicit  dimensioning,  associative  dimensioning,  plus- 
minus  tolerancing,  and  geometric  tolerancing.  This  part  captures  the  semantic  contents  of 
the  dimensions  and  tolerances  as  specified  in  those  standards  and  applies  to  3-D  geometric 
models  of  parts. 

The  major  subdivisions  of  this  part  of  10303  are: 
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shape  aspect  definition:  provides  resources  for  the  representation  of  shapes  to  which 
dimensions  and  tolerances  are  applied; 


• shape  dimension:  provides  resources  for  the  representation  of  size  and  relative  location 
to  meet  the  dimensioning  requirements  found  in  the  engineering  design; 

• shape  tolerance:  provides  resources  for  the  representation  of  limits  within  which  man- 
ufactured shapes  are  permitted  to  vary. 

However,  ISO  10303  does  not  cover  the  following: 

• the  definition  of  the  fundamental  principles,  concepts  and  terminology  of  tolerancing 
and  dimensioning; 

• the  mathematical  definition  of  tolerances  and  datums; 

• the  description  of  dimensioning  or  tolerancing  practices; 

• the  specification  of  dimensioning  inspection  methods; 

• the  synthesis  and  analysis  of  tolerances; 

• the  tolerancing  of  product  characteristics  other  than  shape; 

• the  presentation  of  tolerances  on  engineering  drawings; 

• the  specification  of  surface  finish  or  surface  roughness. 

5.2  ISO  TC213:  Geometrical  Product  Specifications  (GPS) 

ISO  TC  213’s  scope  is  standardization  in  the  field  of  geometrical  product  specifications 
(GPS)  i.e.  macro-  and  micro-geometry  specifications  covering  dimensional  and  geometrical 
tolerancing,  surface  with  properties  and  the  related  verification  principles,  measuring  equip- 
ment and  calibration  requirements  including  the  uncertainty  of  dimensional  and  geometrical 
measurements.  The  standardization  includes  the  basic  layout  and  explanation  of  drawing 
indications  (symbols). 

ISO  TC  213  is  adopting  a bottom-up  approach  for  design.  That  is  starting  from  feature 
level  to  the  product  level.  The  STEP  community  is  adopting  a top-down  approach,  that  is, 
product-assembly-parts-features-conditions.  These  two  approaches  should  logically  meet  at 
some  common  level. 

ISO  5459  [28-30]  is  prepared  by  ISO  TC  213.  ISO  5459  cancels  and  replaces  ISO  5459:1981, 
and  a technical  revision  of  ISO  5449:1981.  It  consists  of  the  following  parts,  under  the  general 
title  Geometrical  product  specifications  (GPS)  - Datums  for  geometrical  tolerancing: 
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Part  1:  Terms  and  definitions 

Part  2:  Datum  and  datum  systems:  Drawing  indication 

Part  3:  Association  methods  for  datums  and  datum  systems  for  the  assessment  of  geomet- 
rical tolerances 

Part  4:  Metrological  establishment  of  datums  and  datum  systems  for  the  assessment  of 
geometrical  tolerances 

The  main  definitions  related  to  5459-1  and  5459-2  are  described  in  ISO/CD  17450:  Ge- 
ometrical Product  specifications  (GPS)-Model  for  geometric  specification  and  verification 
[32].  The  basic  mathematical  idea  is  as  follows:  Given  a set  of  features,  there  are  some 
operations  (partition,  extraction,  collection  etc.)  by  which  we  can  get  some  other  features 
that  belong  to  this  set.  This  way  we  can  build  the  part,  sub-assembly,  assembly,  and  finally 
the  product.  This  is  the  bottom-up  approach  taht  was  mentioned  earlier  in  the  report.  The 
second  part  of  ISO  5459  [29]  describes  the  rules,  explanations  and  manner  in  which  datums, 
and  datum  systems  are  indicated  in  technical  product  specification.  This  part  defines  Single 
datums,  common  datums  and  datum  systems.  It  also  defines  the  writing  and  reading  rules 
of  the  graphical  language  (drawing  indications).  Currently  this  rule  based  system  needs  a 
mathematical  rigor. 

The  WG  14  of  TC  213  is  working  on  ISO/DTS  17450  Geometrical  Product  Specifications 
(GPS)  - Model  for  geometrical  specification  and  verification  [32,33].  This  is  the  basis  on 
which  the  entire  GPS  language  is  built.  This  is  the  most  important  and  fundamental  part. 
It  has  been  decided  to  split  this  standard  into  two  parts 

• ISO/CD  17450-1:  Geometrical  product  specifications  (GPS)-  General  concepts-  Part-1: 
Model  for  geometric  specification  and  verifications; 

• ISO/CD  17450-2:  Geometrical  product  specifications  (GPS)-  General  concepts-  Part-2: 
Operators  and  uncertainties. 

In  this  standard  they  are  classifying  the  various  features  used  in  the  field  of  GPS  and 
defining  the  operators  on  these  features.  For  example,  ideal  feature:  perfect  shape  feature 
defined  by  a type  (geometrical  properties  of  a feature)  and  characteristics  (parameter  of 
one  feature  or  between  two  features  expressed  as  a length  or  an  angle),  non-ideal  feature: 
imperfect  feature  fully  dependent  of  the  skin  model,  and  skin  model:  model  of  the  the  phys- 
ical workpiece  with  its  environment.  The  following  operations  are  defined  in  ISO/CD  17450 
[32,33]. 

Partition:  used  to  identify  bounded  features. 

Extraction:  used  to  identify  finite  number  of  points  from  a feature,  with  specific  rules. 
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Collection:  used  to  consider  some  features  together,  which  together  play  a functional  role. 
It  is  possible  to  build  the  collection  of  ideal  features  or  the  collection  of  non-ideal 
features.  (Composite  feature) 

Construction:  An  operation  called  construction  is  used  to  build  ideal  features  from  other 
features  maintaining  the  constraints. 

WG  14  is  working  on  ” Operator  Principle”  by  which  they  wish  to  study  whether  the  set 
of  features  along  with  these  operators  can  be  given  a good  mathematical  structure.  But  they 
base  their  assumption  that  controlling  geometry  will  imply  controlling  function.  Controlling 
geometry  will  at  the  most  imply  controlling  geometric  part  of  the  function.  Geometry  is  one 
aspect  and  the  other  important  aspect  is  material.  The  second  major  assumption  they  have 
is  that  summation  of  part  function  will  lead  to  product  function,  which  is  not  true  (there 
may  be  some  few  exceptions/chance  cases).  There  are  lot  of  mathematical  and  logical  issued 
still  unresolved. 

The  FAB  architecture  described  in  this  report  can  substantially  help  in  resolving  these 
issues.  Using  the  various  Class  structures  like  Function,  and  Artifact  we  can  clearly  define 
functional  aspects  of  a feature  without  ambiguity.  The  work  of  STEP  part  ISO  10303-47 
[26]  needs  to  be  harmonized  with  the  work  of  ISO  TC  213.  The  FAB  architecture  can  do 
this  very  effectively  by  using  these  class  structures. 

6 Suggestions  for  Future  Research 

In  this  report,  a conceptual  framework  for  representation  of  artifacts,  function,  and  behavior 
have  been  discussed  and  also  issues  related  to  function  to  form  mapping  in  the  conceptual 
design  stages.  However,  the  conceptual  framework  presented  here  needs  further  investigation 
and  analysis  in  respect  to: 

• detailed  function-form-behavior  specification 

• development  of  computational  methodologies 

• further  considerations  in  detailed  tolerance  analysis  and  synthesis 

• design  optimization  considering  various  aspects  affecting  the  product  life  cycle  starting 
from  the  user  specification  to  manufacturing. 

A detailed  specification  and  formal  definition  of  artifacts  in  terms  of  its  parameters  repre- 
senting functional  requirements,  shape,  goal,  constraints,  etc,  would  be  developed  for  putting 
the  concepts  discussed  in  this  report  on  a concrete  footing.  Representation  of  functions  both 
as  abstract  entities  and  subsequent  links  to  physical  objects  with  concrete  shape  and  be- 
havior, needs  more  rigorous  treatment  so  that  formal  qualitative  as  well  as  quantitative 
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reasoning  techniques  can  be  used  to  drive  the  design  process  in  a systematic  manner.  The 
process  of  selection  of  sub-functions  and/or  objects  during  a functional  decomposition  and 
synthesis  process  can  be  based  on  qualitative  reasoning  using  Case-Based  Reasoning  (CBR) 
and  other  methods. 

Further  work  will  be  required  to  carry  out  detailed  development  of  object-oriented  tools 
for  the  integration  of  tolerance  classes  at  the  conceptual  design  stage  and  development  of 
procedures  for  tolerance  synthesis  and  analysis  models. 

The  various  object-oriented  models  mentioned  in  this  report  would  be  used  to  investigate 
the  development  of  an  optimal  product  design  schema  based  on  integration  of  requirements, 
constraints  and  goals  in  all  aspects  of  product  life  cycle  starting  from  user  specification  to 
design,  manufacturability,  after-sales  support  and  recycling,  etc. 

Figure  17  describes  the  overall  future  directions  of  this  research.  The  center  of  it  all  is 
the  Knowledge  Aided  Design.  This  shared  knowledge  base  drives  the  rest  of  the  applications 
and  tools.  This  is  a clear  shift  from  the  present  CAD  systems  which  are  mostly  geometry 
driven  to  more  information  and  knowledge-driven  design  environment.  The  CAD  tools, 
Analysis/Simulation  tools,  Process  planning  tools  (CAPP)  and  NIST  developed  Process 
Specification  Language  (PSL)  [40]  standard  interface,  business  tools  (PDM,  ERP,  Supply 
Chain  Management)  are  all  connected  to  the  central  knowledge  base.  The  various  data 
exchange  standards  (existing,  under  development  and  future  needs)  are  also  illustrated  in 
the  figure.  Apart  from  these  tools  the  Augmented  CAD  tools  which  include  virtual  reality 
(VR)  tools  is  also  included  in  the  future  goal.  Various  researchers  [15, 16,48]  have  developed 
VR  tools  for  assembly /disassembly,  motion  planning  and  simulation  purposes.  The  data 
exchange  mechanisms  and  standards  are  becoming  a major  issue.  In  developing  the  FAB- 
DFT  model  we  have  provided  open  connectivity  to  these  standards. 

7 Conclusions 

A comprehensive  information  model,  FAB  model,  is  described  in  this  report.  This  model  is 
generic  in  the  sense  that  it  can  be  used  for  capturing  other  information  as  well.  The  Unified 
Modeling  Language  (UML)  [8]  model  under  development  will  be  used  as  a proof  of  concept 
and  will  eventually  lead  to  standards.  The  model  can  also  be  used  to  study  the  various 
interoperability  issues  and  manufacturing  system  integration.  In  [53],  a multi-level  approach 
called  the  Design  for  Tolerance  [DFT]  process  was  proposed  which  enables  tolerancing  to  be 
addressed  at  successive  stages  of  design  in  an  incremental  fashion.  The  effective  use  of  the 
FAB  and  DFT  model  for  design  tolerancing,  starting  from  conceptual  stage  of  the  design 
and  continuously  evolving  throughout  the  entire  design  process  to  the  final  detailed  design 
is  also  discussed  in  this  report. 
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Figure  17:  Future  research  directions  - AP:  Application  Protocol 
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Appendix  A 

A Process  Capability  Indices  (PCIs):  A brief  tutorial 

A manufacturing  process  is  characterized  by  the  measurements  of  process  operations  and 
units  output  from  the  process.  In  most  cases  these  measurements  are  assumed  approximately 
normal  so  that  a process  under  statistical  control  may  be  understood  by  first  two  moments 
mean  p and  standard  deviation  a of  the  statistical  distribution. 

A process  specification  usually  consists  of  lower  and  upper  specification  limits  (LSL,  USL ) 
and  a target  T somewhere  between  the  limits  (most  often  T = MSL  = |( LSL  + USL ) which 
is  a symmetric  case).  A capability  index  (random  variable)  is  a unitless  function  of  the  pro- 
cess parameters  (p,a)  and  the  process  specification  (. LSL,T , USL)  designed  as  a tool  to  aid 
in  the  assessment  of  process  performance  [39].  Computation  of  the  capability  index  involves 
finding  the  ratio  of  allowable  process  spread  or  tolerance  and  actual  process  spread  or  vari- 
ation. For  any  particular  process  the  allowable  process  spread  is  a fixed  requirement,  while 
the  actual  process  spread  is  generally  an  unknown  quantity  that  will  have  to  estimated.  The 
two  main  approaches  for  PCIs  are: 

1.  based  on  expected  proportion  of  non-conforming  items  (Cp  and  Cpl t); 

2.  based  primarily  on  loss  function  (Cpm). 

The  first  approach  necessitates  the  assumption  on  the  form  of  the  distribution  while  the 
second  approach  assumes  the  form  of  loss  function.  In  this  report  our  discussion  on  PCIs 
is  limited  to  the  basic  definitions,  properties,  and  usage  as  a quality  tool.  For  a complete 
discussion  about  the  statistical  properties  and  estimation  procedures,  interested  reader  can 
refer  to  [42]  and  the  references  cited  in  that  monograph. 

A.l  The  Basic  PCIs:  Cp  and  Cpk 

Suppose  the  lower  and  upper  specification  limits  for  the  value  of  a measured  characteristic  X 
be  specified  as  LSL  and  USL  respectively.  Values  of  X outside  [LSL, USL]  will  be  termed 
as  non-conforming  and  we  are  interested  in  the  probability 

Pr{LSL  < X < USL}. 

An  indirect  measure  of  potential  ability  (capability)  to  meet  the  requirement  ( LSL  < X < 
USL ) is  the  process  capability  index  Cp: 


67 


Cp  = l 


USL^LSL  Bilateml 
6a 


usf,...i£  Unilateral  with  USL  known 

3(7 

/i~L5L  Unilateral  with  LSL  known 

3(7 


(7) 


where  a denotes  the  standard  deviation  of  X . From  (7)  we  can  clearly  see  the  large  values 
of  Cp  are  desirable.  From  the  definition  of  Cv  we  see 


Specification  width  allowable  process  spread  or  tolerance 
Capability  index  = — — : — — = : : — ; . 

Distribution  width  actual  process  spread  or  variation 

Usually  it  is  assumed  that  manufacturing  process  measurements  follow  normal  distribution, 
that  is,  X ~ N(ii,  a)  and  often  it  is  assumed  that  target  is  the  mid-point  of  the  specification 
limits,  that  is,  T = MSL  — \ [LSL  + USL)  which  is  a symmetric  case. 

The  use  of  factor  6 in  equation  (7)  has  the  following  implication:  if  Cp  = 1 and  the  process 
distribution  is  normal  with  p.  = \[LSL  + USL)  then  the  proportion  of  non-conforming 
products  is  as  small  as  0.27%  which  is  equivalent  to  2700  ppm  (parts  per  million)  and  this 
has  been  a standard  practice  in  statistical  process  control  (SPC). 

If  the  distribution  is  assumed  to  be  normal  and  if  /i  = T = |( LSL  -I-  USL),  then 
the  expected  proportion  of  non-conforming  products  is  given  by  2<F(—  ^),  where  $(.)  is 
the  standard  normal  CDF  (cumulative  density  function)  and  d = (USL~LSL)  is  the  half- 
length  of  the  specification  width.  From  this  we  can  see  that  the  percentage  yield  is  given 
by  100[2<F(3CP)  — 1] . If  Cv  — 1 the  expected  proportion  of  non-conforming  products  is 
0.27%,  which  is  regarded  as  acceptably  small.  For  acceptance  we  should  have  Cv  < c with 
c = 1, 1.33, 1.5,  2.0  corresponding  to  ( USL  — LSL)  = 6 a,  8a,  9 a,  12 a.  It  is  important  to  note 
that  Cp  = 1 does  not  guarantee  that  there  will  be  only  0.27%  of  non-conforming  product.  In 
fact,  all  that  it  does  guarantee  is  that,  with  assumption  of  normality  and  the  relevant  value 
of  a,  there  will  never  be  less  than  0.27%  expected  proportion  of  non-conforming  product. 
It  is  only  when  fi  = T = |(L5L  + USL)  that  the  expected  proportion  of  non-conforming 
product  is  small  as  0.27%.  A good  discussion  of  these  indices  including  estimating  procedures, 
sampling  strategies,  and  application  methodologies  is  given  in  [39]. 


A. 2 Cpk  Index 

It  is  clear  from  the  definition  of  Cp,  it  fails  to  take  into  consideration  the  effect  of  — T\. 
Because  of  inability  of  Cp  to  consider  target,  several  indices  have  been  proposed  that  take 
into  account  the  target  T.  The  Cpk  index  was  introduced  to  give  the  values  of  /z  some 
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influence  on  the  value  of  the  PCI.  Cpk  is  defined  as 

min{USL—n,fi—LSL] 

— 3(7 


d- 


( USL+LSL ) 
3cr 


(since  min(a,  b ) = ^(|a  + b\  — \a  — 6|) 


(8) 


= 1- 


( USL  + LSL ) | 

If4  2 I 


cD 


Since  Cp  = we  have  Cpk  < Cp , with  equality  if  and  only  if  \d  — T — |( LSL+USL ) = 
MSL.  It  seems  that  the  original  way  of  compensating  for  the  lack  of  input  ^ in  the  calculation 
of  Cp,  was  to  use  an  additional  statistic  k (Centering  Index)  defined  as 


k 


\H-MSL\ 

d 

for  any  general  Target  T 


(9) 


This  is  equivalent  to  using  [i  (through  k)  and  a(through  Cp)  separately.  While  it  might 
be  said  that  one  might  as  well  use  /r  and  a themselves  (or  their  estimates),  k and  Cp 
do  have  the  advantage  of  being  dimensionless,  and  related  to  the  specification  limits.  In 
[72],  the  centering  index  k is  denoted  as  Cc.  Johnson  [50]  even  points  out  that  “none  of  the 
indices  adds  any  knowledge  or  understanding  beyond  that  contained  in  the  basic  parameters 
yL,  a , T,  LSL,  USL. 

With  the  combination  into  the  index  Cpk  = (1  — k)Cp,  the  reduction  to  single  index 
indeed  had  to  be  paid  for  by  losing  separate  information  on  location  [yf]  and  dispersion  {a). 
Cpk  and  C*k  are  numerically  equivalent  if  0 < A;  < 1.  Cpk  captures  the  off-targetness.  The 
centering  index  k represents  the  amount  of  tolerance  consumed  by  the  mean  shift. 

Yet  another  way  of  defining  Cpk  is  by  defining  Cpi  and  Cpu  as 


and 


and  define 


^ Id- LSL 

USL-fi 
Upu  ~ 3a 

Cpk  — rmn{Cp/,  Cpu } 


If  the  distribution  is 
product  is 


assumed  to  be  normal,  then  the  expected  proportion  of  non-conforming 


$ 


yi-LSL\ 
a ) 


+ 1 — $ 


USL- 


and  we  have 


100[$(3Cpfc)  - 1]  < %Yield  < 100[$(3Cpfc)]. 
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A. 3 Cpm  Index 


The  Cp  and  CPk  indices  are  appropriate  measures  of  progress  for  quality  improvement  in 
which  ‘reduction  of  variability’  is  the  main  criterion  and  process  yield  is  the  primary  measure 
of  success.  Taguchi  [78,  79]  has  suggested  a different  approach  to  quality  improvement  in 
which  ‘reduction  of  variation  from  the  target  value’  is  the  guiding  principle.  The  index  CPm 
is  defined  as 


a 


pm 


USL— LSL 

f>yJa2+(p-T)2 


(10) 

d 

3^/a2+(n-T)2 


The  cost  of  a characteristic  X missing  the  target  T is  often  assumed  to  be  well  approxi- 
mated by  the  symmetric  squared  error  loss  function 


loss{X ) = c(A-T)2 

where  c is  some  positive  constant.  Since  E[(X  — T)2]  = E(X2)  + T2  - 2 pT  -I - p2  — p2  = 
a2  + (p  — T)2  = E[loss(X)],  the  Cpm  index  can  be  written  as 


^ _ USL  - LSL 

Is  pm  — / 

6 yJE[loss(X)\ 

The  relationship  of  CpTn  to  squared  error  loss  is  discussed  in  [50].  Chan  et  al  [11]  proposed 
the  index  Cpm  with  the  assumption  p = T.  Boyles  [10]  extended  the  definition  without  the 
restrictive  assumption  p = T and  explains  how  a given  value  of  Cpm  places  an  upper-bound 
on  |p  — T|. 

Since  a2  + (p  — T)2  > a2  we  have  Cpm  < Cp  and  Cpm  = Cp  <=>  p = T.  The  index  Cpm 
as  defined  is  not  a satisfactory  measure  of  process  capability  unless  the  target  value  is  equal 
to  the  mid-point  of  the  specification  limit  [42].  A modified  index  Cpm  is  defined  to  overcome 
this  problem. 

^ min  {USL-T,T-LSL} 

3yV  + („  - Tf 

If  r = MSL  then  Cpm  = C;m. 


A. 4 Comparison  of  PCIs 

If  p = T,  then  Cp  = CPk  — Cpm ■ If  the  process  distribution  ~ A(p,  a)  (and  in  statistical 
control)  then  a value  of  1 for  all  the  PCIs  means  that  the  expected  proportion  of  product 
within  the  specification  limits  is  > 99.7%.  Negative  values  are  not  permitted  for  PCIs. 

As  already  explained,  the  main  function  of  Cpk  is  to  indicate  the  degree  to  which  the 
process  is  within  specification  limits.  For  LSL  < p < USL , CPk  — > oo  as  a — » 0,  but  large 
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values  of  Cpm  do  not  provide  information  about  \fi  — T\.  The  index  Cpm  — » oo  if  /i  — > T as 
well  as  a — * 0 and  Cpm  — » 0 as  \fi  — T\  — » oo.  The  index  Cpm  is  bounded  above  as 


^ ^ USL  - LSL 

C pm  < 71  77 

6|/i  — 1 I 


(11) 


It  is  clear  from  (11),  when  T = MSL , a Cpm  value  of  1 implies  that  the  process  mean  /i  lies 
within  the  middle  third  of  the  specification  width. 

In  [69]  a unifying  approach  is  discussed  that  ties  the  various  indices  together.  Through 
the  use  of  a weight  function  the  relationship  that  exists  among  Cp,  CPk , Cpm  and  Cpk  is 
discussed  [69].  The  index  Cpw  is  defined  as 


C 

^ pw  — 


USL  - LSL 
6^/cr2  + w(fi  — T): 


(12) 


where  w represents  the  weight  function,  Cp0  = Cp  and  Cpi  = Cpm.  For  complete  details  refer 
[69].  There  are  other  additional  weight  functions,  for  example, 


w = 


r 0 < r 
0 otherwise 


(13) 


- 0 < r 

T 

0 otherwise 


(14) 


The  weight  function  (13)  would  typically  be  used  in  those  situations  where  minor  depar- 
tures from  the  target  are  not  of  interest,  but  major  departures  are  considered  critical.  The 
weight  function  (14)  applies  greater  weight  to  minor  departures  from  the  target,  while  only 
a marginal  decrease  in  capability  results  once  the  process  deviates  from  the  target  by  more 
than  la.  Using  different  weight  functions  allows  one  to  customize  the  PCI  to  the  process  of 
interest,  thereby  allowing  different  shaped  loss  functions  (as  in  the  case  of  Cpm). 


A.  5 Example 

Suppose  we  have  three  different  lots  of  fasteners  (not  measured)  and  we  are  interested  in 
their  values  (say,  length).  We  can  compute  the  sample  mean  x — where  n is  the 

sample  size.  We  assume,  for  the  sake  of  simplifying  our  discussion,  that  sample  mean  and 
variance  is  equal  to  population  mean  and  variance.  This  implies  that  our  sample  is  a perfect 
representation  of  the  population,  which  we  assume  as  a Gaussian  population. 

Let  the  specification  limits  be  specified  as  LSL  = 114  mm,  USL  = 126  mm  and  the 
target  T — 120  mm.  Let  fix  = 120 ,/i2  — 123,  jU3  = 124  and  o\  — 2,a2  = 1 , cr3  = 0.6666 
be  the  means  and  standard  deviations  of  three  populations  respectively.  From  the  theory 
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of  the  normal  distribution,  we  can  compute  the  following  probabilities  (Lot#l  is  taken  as 
example). 


[i  ± 1(7 

68.26% 

Prjfastener  £ [118,122]}  = 0.6826 

[i±2a 

95.46% 

Prjfastener  £ [116,124]}  = 0.9546 

fi  ± 3a 

99.73% 

Pr{fastener  £ [114,126]}  = 0.9973 

fi±6cr 

99.9999998% 

Prjfastener  £ [108,132]}  = .999999998 

From  the  table  we  can  see  that  99.73%  of  the  fasteners  lie  in  [114,126]  mm.  The  expected 
proportion  of  non-conforming  products  is  0.27%  or  the  defect  rate  is  2700  ppm  (parts  per 
million).  [2700  ppm  is  computed  as  (1  — 0.9973)  x 1000000].  The  defect  rate  for  Lot#l 
would  be  0.002  ppm  if  we  have  Let  us  now  compute  the  PCIs  Cp,  Cpk  and  Cprn. 


Lot#  1 

Lot#  2 

Lot#  3 

k 

0 

|123:120'  = o.5 

D 

1 124 — 120|  _ 2 

6 ~ 3 

cp 

1 

2 

3 

Cpk 

1 

1 

1 

0 

3 

6 

C — 1 

1 

0.63 

0.49 

Note  that  Cpm  clearly  distinguishes  the  three  different  lots.  The  index  Cp  is  increasing 
from  Lot#l  to  Lot#3  while  Cpm  is  decreasing  from  Lot#l  to  Lot#3.  It  is  important  to  note 
that  the  mean  is  drifting  from  the  target  for  lots  two  and  three.  The  index  Cpk  is  equal  to 
one  for  all  the  lots  and  does  not  capture  process  centering. 

A. 6 Generalization  of  PCIs 

The  previous  discussion  about  PCIs  assume  the  normal  distribution  for  the  process  (or  for 
the  measured  characteristic  X),  which  is  a strong  assumption.  We  would  like  to  drop  this 
normality  assumption.  This  may  be  necessary,  because  perfect  normality  is  seldom  realized 
in  practice,  and  even  in  theory  the  normal  distribution  might  not  be  expected  for  certain 
situations  (for  example,  perpendicularity  tolerance). 

The  study  of  non-normal  distribution  for  X falls  into  two  main  parts.  The  first  is  the 
investigation  of  the  properties  of  PCIs  and  their  statistical  estimates  when  the  distribution 
has  specific  non-normal  forms.  The  second,  more  difficult,  is  development  of  methods  of 
allowing  for  non-normality  and  consideration  of  use  of  new  PCIs  specifically  designed  to  be 
robust  (that  is,  not  too  sensitive)  to  non-normality.  There  are  several  methods  suggested  in 
the  literature,  here  we  discuss  briefly  about  (1)  Clements’  Method  and  (2)  Johnson-Kotz- 
Pearn  Method.  For  complete  details  refer  [42]. 
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A. 6.1  Clements’  Method 


Clements  [14]  proposed  a method  of  construction  based  on  the  assumption  that  the  process 
distribution  can  adequately  be  represented  by  Pearson  distribution.  His  suggestion  was  to 
replace  the  multiplier  ‘6’  in  the  denominator  of  Cp  by  number  9 such  that 

Pr  j/i  — -6o  < X < p + —do | = 0.0027. 

For  given  values  of  skewness  (v^i)  and  kurtosis  (/52)  coefficients,  compute  9 such 
the  above  probability  is  true.  For  example,  if  y/]T]  = 1 and  f32  = 5,  we  have  6 = 6.572 
The  index  Cp  would  be  calculated  as  y^-.  For  normal  distribution  we  have  9 = 6 so 
Cp  = For  complete  discussion  about  this  methods  refer  [14]. 

A. 6. 2 Johnson-Kotz-Pearn  (JKP)  Method 

Clements’  method  requires  the  evaluation  (estimation)  of  and  (32  which  may  be  difficult 
to  obtain.  We  need  rather  large  samples  for  accurate  estimation  these  quantities.  This 
method  defines  a new  PCI  Cp^  as  [42] 

^ _ USL  - LSL  _ d 

Cm  = — — = (1^) 

where  9 is  chosen  so  that  the  ‘capability’  - namely  the  proportion  of  conforming  items,  with 
optimum  choice  of  p is  not  greatly  affected  by  the  shape  of  the  distribution.  Note  that 
Cp  = Cp(e= 6)-  A value  of  5.15  is  recommended  for  9 [42],  so  that  Cp  = Cp^=5. 15)  = 0 ft.oa  and 

r _ d-\n-MSL\ 

L/Pk  ~ 2. 5750- 

Note  that  robustness  is  attained  at  the  cost  of  abandoning  the  0.27%  level  for  expected 
proportion  of  non-conforming  items  (2700  ppm).  Both  methods  rely  on  the  assumption  that 
the  population  distribution  is  uni-modal  shape  close  to  a Pearson  distribution  for  Clements 
method,  and  more  restrictively,  close  to  a Gamma  distribution  for  JKP  method. 

Remark  A.l  The  natural  variation  (shifts  and  drifts,  for  example,  due  to  tool  wear , 
temperature  etc)  can  cause  the  mean  to  deviate.  Motorola  suggests  the  mean  shift  could  be  of 
the  order  of  1.5o  units  from  its  target  [24] ■ For  a product  to  be  Motorola  six  sigma  quality, 
that  is  d — 6 a,  it  has  to  have  Cp  > 2,  CPk  > 1.5  with  a defect  rate  of  3.4ppm.  Section  2.3.1 
(page:  19)  discuss  this  in  detail. 

The  Motorola  convention  is  to  use  a one-sided  mean  shift  of  1.5cr.  The  one-sided  mean 
shift  is  perhaps  motivated  by  common  physical  phenomena  such  as  tool  wear.  A shift  o/1.5a 
is  motivated  by  earlier  work  by  Bender  [6].  Also,  it  is  assumed  that  the  process  standard 
deviation  is  invariant  (Appendix  B,page:  76). 


that 

[14]. 

that 
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If  Cp  — 2 and  CPk  = 1.5  (mean  shift  consumes  25  percent  of  the  tolerance  range),  the 
probability  of  conformance  can  be  shown  to  be  0.9999966,  which  is  equivalent  to  3.4  ppm. 
Thus  Cp  > 2 and  Cpk  >1.5  imply  six  sigma  quality,  assuming  a 1.5 a one  sided  mean  shift. 

The  table  A. 6. 2 summarize  the  various  capability  indices  we  have  discussed. 


Remark  A. 2 The  ISO  effort  on  standards  for  statistical  tolerancing  [72,  73]  describes 
the  acceptable  statistical  distributions  of  the  population  of  features  using  Cp,Cpk  and  Cpm- 
At  a conceptual  level  statistical  tolerancing  simply  defines  a class  of  acceptable  distributions 
of  the  population  of  features  but  for  practical  applications  the  methods  to  define  and  compute 
these  classes  should  be  provided.  In  [73]  three  different  methods  of  defining  the  classes  are 
explained:  (1)  characterizing  the  distributions  by  their  first  two  moments;  (2)  Cp  — Cpk  plot; 
and  (3)  CDF  based.  Refer  Section  2.5,  page: 25. 

Remark  A. 3 All  process  capability  indices  of  this  from  are  determined  under  the  as- 
sumption that  (1)  the  process  is  in  statistical  control,  (2)  the  target  is  the  mid-point  of  the 
specification  limits,  and  (3)  the  process  measurements  are  normally  distributed.  If  the  pro- 
cess is  centered  and  the  target  is  the  mid-point  of  the  specification  limits,  then  the  indices  are 
equivalent.  If  the  last  two  assumptions  are  relaxed  (that  is  non-symmetric  and  non-normal) 
care  should  be  taken  in  interpreting  the  values  of  the  indices.  But  as  long  as  the  distributions 
are  determined  by  their  first  two  moments  these  indices  will  serve  as  an  useful  tool  to  assess 
the  process  performance. 
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LSL,  USL 
MSL 


Lower  and  Upper  Specification  Limit 
MSL  = 

Mean  and  standard  deviation  of  the  process  distribution 


T 


Target  specification 


j _ USL— LSL 
u — 2 


Index  Definition 


Properties 


Cn 


USL-LSL 

6 a 


Bilateral 


0SJff  ^ Unilateral  USL  known 
|U~3^sl  Unilateral  LSL  known 


X ~ N(iL.cr),  n = T = \{LSL  + USL ) 
Cp  = 1 =>  expected  proportion 
of  non  — conforming  products  < 0.27% 
Percentage  Yield  = 100[2$(3Cp)  - 1] 

Cp  ^ max{Cpyt,  CPm} 


a 


pk 


mm 


{ USL—p,p—LSL } 


3 a 


, I ( USL+LSL ) 

d-\p-  2 

3a 


I ( USL+LSL ) I 

= [)Cp 


Cpk  Si  Cp  and  Cpk  - Cp  <==$■  /i  = T = MSL 
Cpk  — 0 if  RHS  is  — ve 
IfX  ~ N(fjL,  cr),fj,  = T = MSL 
100[2$(3CP^)  - 1]  < %Yield  < 100[<L(3Cpfc)] 


k 


\p-MSL\ 


I H-T\ 
~ d 


Centering  Index 


for  any  general  Target  T 


C. 


pk 


(1  - k)Cp 


0 < k < 1 =>  Cpk  — C. 


''pk 


c 


pm 


USL-LSL 

6y/a2+(p-T)2 

d 


r < r 

Ky'pm  — ^ p 


3y/a2+(»-T)2 

USL-LSL 

6y/E[loss(X )] 


Cpm  — Cp 

Cpm  as  defined  is  not  a statistical 


[i  = T = MSL 
; not  a statistica 
measure  of  process  capability  unless  T — MSL 


Table  2:  PCIs  Summary 
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Appendix  B 

B The  Defect  Rate  of  3.4  ppm  and  Motorola  Conven- 
tion 

A part  or  an  item  is  classified  as  defective  (non-conforming)  if  the  desired  measurement 
(quality  characteristic),  denoted  by  A,  is  outside  the  customer’s  upper  specification  limit 
(USL)  or  lower  specification  limit  (LSL).  In  addition  to  specifying  the  USL  and  LSL,  a 
customer  would  also  specify  a target  value,  T,  which  is  typically  the  midpoint  of  specification 
width.  For  example,  let  us  say  that  the  customer’s  specification  on  the  length  of  fasteners  is 
120  ± 12  mm  ( USL  — 132,  LSL  = 108,  T = 120  and  the  half-tolerance  width  is  12mm).  The 
process  that  produce  these  fasteners  is  such  that  the  lengths  are  normally  distributed.  We 
assume,  for  the  sake  of  simplifying  our  discussion,  that  sample  mean  and  variance  is  equal 
to  population  mean  and  variance.  This  implies  that  our  sample  is  perfect  representation  of 
the  population,  which  we  assumed  it  as  normal  population. 

Assuming  that  the  process  is  centered,  the  number  of  defects  per  million  for  the  different 
quality  levels  is  given  in  Table  3. 


Quality  Level  (a) 

Defect  Rate  (ppm) 

1 

317310 

2 

45500 

3 

2700 

3.5 

465 

4 

63 

4.25 

22 

4.5 

6.8 

4.75 

2.0 

5 

0.6004 

5.25 

0.1636 

5.5 

0.0424 

5.75 

0.01524 

6 

0.002 

Table  3:  Quality  levels  and  defect  rates 

As  table  3 shows,  there  are  2700  defects  per  million  in  the  3 a quality  level  and  only  63 
defects  per  million  parts  in  the  4a  quality  level  and  only  2 defects  per  billion  parts  in  the  60- 
quality  level.  But  then  how  do  we  get  the  defect  rate  of  3.4  ppm  for  6o-  quality  level?  We 
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will  now  show  the  actual  computation  for  the  6a  quality  level.  The  figures  in  the  above  table 
are  arrived  at  using  the  normal  probability  tables  under  the  assumption  that  the  process  is 
centered.  Refer  figure  18  to  understand  the  typical  areas  under  the  normal  curve.  The  area 
under  the  pobabilty  curve  for  any  general  distribution  is  given  by  Tchebychev’s  inequality. 
If  the  process  mean  is  equal  to  T = |( USL  4-  LSL ),  the  process  is  centered;  otherwise  it  is 


— 8cr  — 7 cr  — 6cr  -5  a — -4cr  — 3cr  — 2c?  — Icy  ju.  +la  -1-2  cr  -t-3cr  -4-4-ct  -i-5ct  -l-6cr  — 8cr  —Icy 

Figure  18:  Areas  under  the  normal  curve 

off-centered.  Because  of  typical  shifts  and  drifts  in  the  process,  the  process  is  off-centered 
and  the  process  mean  shifts  by  as  much  as  1.5a  units  [22].  A shift  of  1.5a  is  motivated  by 
earlier  -work  by  Bender  [6]  and  Gilson  [19].  More  about  this  later. 

The  quality  level  can  be  expressed  as  ma,  where  m — ^ S-L~a . Thus  if  a = 4 (Cp  — 1), 
the  process  has  3a  quality  level;  if  a = 3 (Cp  = 1.33),  the  process  has  4a  quality  level;  and 
if  a = 2 (Cp  = 2),  the  process  has  6a  quality  level. 

Motorola  has  established  a convention  of  distributing  the  total  proportion  of  the  non- 
conformity to  the  right  side  of  the  normal  curve.  This  one-sided  mean  shift  is  perhaps 
motivated  by  common  physical  phenomena  such  as  tool  wear  taht  becomes  increasingly 
worse.  Also,  it  is  assumed  that  the  process  standard  deviation  is  invariant. 

The  probability  calculations  with  respect  to  actual  are  (Figure  19) 

1 — Pr  {conformance} 

= 1 - Pr{XAe[LSL,USL}} 

= 1 - Pr{  108  < < 132} 

= 1 - Pr{123  - 15  < AT  < 123  + 9} 
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= 1 — Pr{123  — 7.5a  < Xa  < 123  + 4.5a} 
- 1 — Pr{— oo  < Xa  < p.  + 4.5a} 

= 1 - 0.9999966 
= 3.4ppm  defect  rate 


As  said  earlier,  if  Cp  = 2 and  = 1.5  (mean  shift  consumes  25  percent  of  the  tolerance 
range),  the  probability  of  conformance  can  be  shown  to  be  0.9999966,  which  is  equivalent  to 
3.4  ppm  defect  rate.  Thus  Cv  > 2 and  Cpk  > 1.5  imply  six  sigma  quality,  assuming  a 1.5a 
one  sided  mean  shift. 

In  the  example,  this  means  that,  if  the  process  producing  the  fasteners  has  a six-sigma 
quality  level  (a  = 2)  and  if  the  process  mean  is  off-centered  by  as  much  as  1.5a  units 
(117  < \i  < 123),  the  maximum  number  of  defects  is  3.4  per  million. 

B.l  The  1.5(j  Mean  Shift 

As  explained  in  the  previous  section  the  six-sigma  quality  implicitly  assumes  that  the  process 
mean  can  be  off-centered  by  as  much  as  1.5a.  It  is  very  difficult  to  accept  that  any  process 
would  drift  as  much  as  1.5a  without  being  noticed  or  corrected. 

According  to  Bender  [6],  the  1.5a  mean  shift  is  based  on  “probability,  approximation  and 
experience”  considerations  and  sometimes  referred  as  Benderizing.  Unfortunately,  Bender 
and  Gilson  introduced  the  1.5a  factor  in  the  different  context  of  tolerancing  assemblies  of 
components.  Evans  [17],  in  his  series  of  articles,  discusses  this  concept  with  a classic  example 
of  stack  of  n nominally  identical  disks.  For  illustrative  purposes,  consider  n = 10  and  the 

height  of  the  stack  as  the  response  (assembly  response  function  Y = X\  -I-  A2H \-Xn  [53]). 

Let  us  assume  that  the  stack  height  must  be  1.25  ± 0.01  in.  The  problem  with  statistically 
tolerancing  the  components  of  an  assembly  (individual  disks)  deals  with  setting  limits  on 
the  allowed  excursion  of  the  components’  mean  and  range  (using  the  standard  deviation).  If 
a2  is  the  variance  of  the  component  z,  the  variance  of  assembly  is  given  by 

a2  = Var (assembly)  = a 2 — na2  if  a*  = ac;  Vz  (15) 

Traditionally,  the  half  width  of  the  specified  tolerance  is  set  to  three  times  the  standard 
deviation  of  the  assembly  in  determining  the  standard  deviation  of  each  component.  In  this 
example: 

0.01  = 3aa  = 3\/l0ac 

Thus,  the  process  would  be  designed  so  that: 

0.01 
ac“  Sv/10 
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Figure  19:  Distribution  with  1.5cr  mean  shift 


Bender  suggests,  however,  that  due  to  potential  shifts  in  the  means  of  the  disks,  it  is  practical 
to  use  1.5aa  as  the  effective  standard  deviation  of  the  assembly.  Thus,  he  recommends  using 
aa  = 1.5  ^2  of  in  (1).  Harry  [22]  seems  to  interpret  the  factor  of  1.5  as  the  allowance  for 
the  shifts  in  the  mean  of  a single  component,  due  to  its  being  manufactured  in  different 
lots.  As  suggested  by  Bender  and  Gilson,  it  might  be  reasonable  to  inflate  the  estimator  of 
the  standard  deviation  of  the  assembly  to  allow  for  the  shifts  and  drifts  in  the  mean  of  the 
individual  components.  However,  we  must  justify  the  1.5cr  shift  in  the  process  mean  of  the 
individual  components. 

In  many  situations,  adjusting  the  process  to  move  the  process  mean  closer  to  the  target 
value  is  relatively  easier  and  cost-effective  than  improving  the  process  to  reduce  the  variance. 
Thus,  if  the  goal  is  to  reduce  the  number  of  defects,  it  does  not  make  sense  to  improve  the 
process  to  six-sigma  levels  and  not  center  the  process. 

It  is  interesting  to  note  that  the  desired  quality  levels  (expressed  by  the  number  of 
defects  in  ppm)  might  be  achieved  through  several  combinations  of  off-centering  and  process 
standard  deviation.  For  example,  a quality  level  of  3.4  defects  per  million  parts  can  be 
achieved  in  atleast  three  different  ways:  0.5a  off-centering  with  5a  quality;  la  off-centering 
with  5.5a  quality;  and  1.5a  off-centering  with  6a  quality.  How  to  achieve  a specified  quality 
level  or  a given  number  of  defects  per  million  depends  on  the  costs  associated  with  adjusting 
the  process  mean  versus  reducing  the  process  variance.  If  the  process  mean  can  be  centered, 
monitored,  and  maintained  at  the  target  value,  a 4a  quality  level  results  in  only  63  defects 
per  million,  and  even  a 3.5a  program  results  in  just  465  defects  per  million  parts.  Table  4 
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Quality 

level 

Mean  Shift  (in  a units) 

0 

0.25 

0.5 

0.75 

1.0 

1.25 

1.5 

1.75 

2 

3.0a 

3.5a 

4.0a 

4.5a 

5.0a 

5.5a 

6.0a 

2700 

465 

63 

6.8 

0.57 

0.034 

0.002 

3577 

666 

99 

12.8 

1.02 

.1056 

0.0063 

6440 

1382 

236 

32 

3.4 

.71 

0.019 

12288 

3011 

665 

88.5 

11 

1.02 

.1 

22832 

6433 

1350 

233 

32 

3.4 

.39 

40111 

12201 

3000 

577 

88.5 

10.7 

1 

66803 

22800 

6200 

1350 

233 

32 

3.4 

105601 

40100 

12200 

3000 

577 

88.5 

11 

158700 

66800 

22800 

6200 

1300 

233 

32 

Table  4:  The  number  of  defectives  (ppm)  for  specified  process  mean  shift  and  Quality  levels 


gives  the  number  of  nonconforming  parts  per  million  for  different  quality  levels  and  different 
values  of  mean  shift.  If  the  process  centering  cannot  be  monitored  or  effectively  controlled,  a 
little  latitude  would  be  provided  on  each  side  of  the  specification  to  have  some  guard  against 
process  shifts. 
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